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ABSTRACT 

The  Mtissbauer  effect    has   been  vised  to  investigate  electronic 
shielding  by  closed  electron  shells   in  salts   of  trivalent  thulium,    by 

measuring  the  temperature  dependence  of  the  nuclear  quadrupole 

169 
splitting  of  the   8.  42  keV  gamma  transition  in  Tm         .      The  nuclear 

quadrupole  interaction  was   studied  for   Tm         ions   in  thulium  ethyl 

sulfate,    thulium  oxide  and  thulium  trifluoride  within  a  temperature 

range  from  9.6    K  to  1970    K.     The  interpretation  of  the  experimental 

data  in  terms   of  the  contributions   of  distorted  closed  electron  shells 

to  the  quadrupole  interaction  yields  values   for  electronic   shielding 

factors.      The  results   lead  to  amounts   of   10%  or  less  for  the  atomic 

Sternheimer  factor  R^.      The  experiments  also  reveal  substantml 

shielding  of  the  4f  electrons  from  the  crystal  electric  field,    ex  - 

pressed  by  the  shielding  factor  cr    .      Values   of  250  and   128  are  ob- 

tained  for  the  ratio  ( 1  —  y     )/(l-(r,)  for  thulium  ethyl  sulfate  and  thu- 

'co  2  J 

lium  oxide  respectively,    where  y       is  the  lattice  Sternheimer  factor. 
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I.     INTRODUCTION 

The  technique  of  recoilless   nuclear  resonance  absorption  of 
gamma  radiation,    the  so-called  MHssbauer  effect  (1),    has  been  em- 
ployed in  numerous  experiments   in   recent  years  (2),    (3).      By  binding 
a  radioactive  nucleus  in  a  crystal  lattice  the  emitted  gamma  radiation 
will,    under  certain  conditions  (2),    have  essentially  the  natural  line 
width  as  determined  by  the  Heisenberg  uncertainty  relation  and  an 
energy  exactly  equal  to  the  excitation  energy  of  the  nucleus.     If  a 
nucleus  of  the  same  isotope  which  is  in  its  ground  state  is  also  bound 
in  a  lattice,    there  is   a  large  probability  for  nuclear   resonance  ab- 
sorption of  the  gamma  radiation.    A  distinct  advantage  of  this  tech- 
nique is  the  inherently  high  energy  resolution  that  is  available.      For 
example  in  the  experiments  to  be  described  here,    the  resolution  is 
one  part  in   10       .      Energy   resolutions  of  this  order  make  it  possible 
to  study  nuclear  properties  as  well  as  solid  state  effects  in  the  crys- 
tals that  are  used  to  bind  the  nuclei.      We  make  use  of  the  MHssbauer 
effect    here  to  study  the  nuclear  hyperfine  interactions  in  salts  of 
rare  earths,    specifically  thulium  salts. 

Measurements  of  the  nuclear  quadrupole  interaction  in  salts  of  the 
rare  earth  elements  yield  information  on  the  quadrupole  moments  of 
the   relevant  nuclear  states  and  on  the  electric  field  gradients  which 


(1)  R.  L.    MHssbauer,    Z.    Physik    151,    1  24  ( 1958) ;  Naturwissen- 
schaften_45,    538  (1958);   Z.    Naturforsch.     14a,    211(1959) 

(2)  See  for  instance  H.    Frauenfelder ,    The  MHssbauer  Effect, 
(W.A.    Benjamin  Inc.  ,    New  York,     1962) 

(3)  The  Proceedings  of  the  Third  International  Conference  on  the 
MHssbauer  Effect  appear     in  Rev.  Mod.  Phys .   3_6,  333-504  (1964) 
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exist  in  the  salts  at  the  nuclear  sites.      The  extraction  of  the  compo- 
nents of  the   electric  field  gradient  tensor  from  such  measurements 
is   rather   straightforward  if  the  values  of  the  nuclear  quadrupole  mo- 
ments have  been  obtained  by  other  methods   such  as   Coulomb  excita- 
tion techniques.      On  the  other  hand  the  determination  of  nuclear  mo- 
ments of  rare  earth  nuclei  by  measurements  of  the  nuclear  quadru- 
pole interaction  is   rather  involved  since  this  requires  a  calculation 
of  the  components   of  the  electric  field  gradient  tensor  at  the  nuclear 
sites.      A  calculation  of  the  electric  field  gradients  for  salts  of  the 
rare  earths  can  be  performed  at  present  only  with  limited  accuracy. 
Uncertainties  in  excess  of  30%  are  typical.      It  therefore  appears 
that  measurements  of  the  nuclear  quadrupole  interaction  in  solids  of 
the  rare  earths  are  at  present  of  more  importance  for  studies  of  the 
sources  of  the  electric  field  gradients  than  for  determination  of 
nuclear  quadrupole  moments. 

The  electric  field  gradient  at  the  nuclear  site  of  a  certain  ion 
originates  from  a  number  of  different  sources.      Major  sources  are 
distortions  of  the  electronic  shells  of  the  ion.     These  distortions 
result  from  the  interactions   of  the  electrons  of  the  ion  with  the  crystal 
electric  field  (CEF)  produced  by  the   surrounding  ions  in  the  lattice, 
provided  the  arrangement  of  the   surrounding  ions   reflects  a  point 
symmetry  lower  than  cubic.      The  field  gradient  at  the  nuclear   site 
results   not  only  from  the  distorted  partially  filled  4f  electron  shell 
of  the  rare  earth  ion,    but  also  from  distorted  closed  electron  shells. 
These  distortions   of  the  closed  electron  shells  of  the  ion  constitute  a 
major  source  of  uncertainty  in  calculations  of  the  electric  field 
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gradient  at  the  nuclear  site.      The  deviations  of  the  closed  shells  from 
spherical  symmetry  (electric   multipole  polarization)  usually  lead  to 
substantial  reduction  or  enhancement  (shielding  or  antishielding)  of 
the  electric  field  gradient  at  the  nuclear  site.     Sternheimer  (4),  (5) 
was  first  to  emphasize  the  importance  of  magnetic  dipole  and  electric 
quadrupole  polarizations  of  closed  shells  and  pioneered  in  calculating 
the  contributions  of  closed  shell  deformations  to  the  nuclear  hyper- 
fine  interactions. 

The  nuclear  quadrupole  interaction  depends  strongly  on  the  elec- 
tronic  state  of  the  ion.      The  electronic  states  which  arise  when  a 
rare  earth  ion  is  incorporated  in  a  crystal  lattice  are  basically 
caused  by  the  interaction  of  the  CEF  and  the  electrons  in  the  partially 
filled  4f  electron  shell,    but  the   splittings   of  these  electronic   levels 
are  also  strongly  influenced  by  distortions   of  the  closed  electron 
shells  (6),  (7),  (8).     In  order  to  account  for  the  modification  of  the 
CEF    splitting  which  results  from  electronic   shielding,    one  has  to 
consider  the  quadrupole  moment  as  well  as   higher  multipole  mo- 
ments induced  in  the  closed  shells. 


(4)  R.M.    Sternheimer,    Phys .    Rev.    80,    102(1950);   105  ,    158 
(1957);     R.M.    Sternheimer  and  H.  M.    Foley,    ibid.    102,    731 
(1956);     H.M.    Foley,    R.M.    Sternheimer  and  D.    Tyko,    ibid. 
9_3.    734  (1954) 

(5)  R.M.    Sternheimer,    Phys.    Rev.    84,    244(1951);   _9_5  ,    736(1954) 

(6)  D.T.    Edmonds,    Phys.    Rev.    Letters  J_0,     129(1963) 

(7)  R.G.    Barnes,    E.    Kankeleit,    R.  L.    Mbssbauer  and  J.M. 
Poindexter,    Phys.    Rev.    Letters   11,    253(1963) 

(8)  J.    Blok  and  D.  A.    Shirley  (private  communication) 
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Rare  earth  ions   exhibit  CEF  splittings  which  are  usually  very- 
much  smaller  than  similar   splittings   observed  for  ions   of  the  iron 
transitions  elements.      In  the  iron  transition  series,    the  partially 
filled  3d  electron  shell  is  fully  exposed  to  the  CEF  produced  by 
surrounding  ions,    resulting  in  large  CEF  level  splittings.      The   rela- 
tively small  CEF  level  splittings  observed  for  rare  earth  ions, 
which  typically  are  of  the  order  of  a  few  hundred  cm       ,    probably 
arise  because  of  large  shielding  effects   resulting  from  the   5s    p 
electronic   shells  which  surround  the  partially  filled  4f  shell. 

Present  theoretical  predictions   of  the  influence  of  electronic 
shielding  upon  the  CEF  level  splitting  of  rare  earth  electronic   levels 
diverge.      Burns  (9)  concluded  that  electronic   shielding  in  the   rare 
earth  ions  is  of  little  importance  and  that  the  difference  between  the 
CEF  level  splittings  in  the  iron  series  and  those  in  the  rare  earth 
series   cannot  be  attributed  to  electronic   shielding  of  the  4f  electrons 
from  the  CEF  by  outer  closed  electron  shells.      In  contrast,   Lenander 
and  Wong  (10),    Ray  (11)  and  Watson  and  Freeman  (12)  conclude  that 
electronic  shielding  plays  a  significant  role  in  rare  earth  CEF  level 
splittings . 

Quantitative  estimates  of  actual  shielding  effects  are  hampered 
by  the  lack  of  sufficiently  accurate  atomic  wave  functions  for  rare 


(9)  G.    Burns,    Phys.    Rev.     LZJ^,    2121   (1962) 

(10)  C.J.    Lenander  and  E.  Y.    Wong,    J .  Chem.  Phys .   3_8,    2750(1963) 

(11)  D.K.    Ray,    Proc.    Phys.    Soc.    82_,    47(1963) 

(12)  R.E.    Watson  and  A.  J.    Freeman,    Phys.    Rev.    133,    A1571  (  1964) 
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earth  ions.      Inadequate  knowledge  of  the  contributions  of  the  core 
electrons   is   a  primary  source  of  uncertainty  in  our  present  under- 
standing of  hyperfine  interactions  in  rare  earth  (as  well  as  in  other) 
elements.      Direct  measurements  of  the  influence  of  electronic 
shielding  upon  the  nuclear  hyperfine  interactions  and  upon  the  CEF 
splittings  of  electronic  levels  therefore  are  highly  desirable. 

This  paper  demonstrates  the  use  of  the  technique  of  recoilless 
nuclear  resonance  absorption  of  gamma  radiation  as  a  means  to  ob- 
tain information  on  electronic  shielding  effects  in  rare  earth  isotopes. 
The  procedure  introduced  here  consists  of  combining  measurements 
of  the  temperature  dependent  nuclear  quadrupole  interaction  (per- 
formed by  using  the  technique  of  recoilless   resonance  absorption) 
with  measurements   of  the  CEF  level  splittings  (performed  by  using 
optical  techniques).     Specifically  we  report  on  determinations  of  the 
relevant  electronic  shielding  factors  for  trivalent  thulium  based 

upon  our   gamma-absorption  measurements   of  the  nuclear  quadrupole 

1  69 
interaction  of  Tm  in  thulium  ethyl  sulfate  (13)  and  thulium  oxide 

and  on  optical  measurements  of  CEF  levels  by  Wong  and  Richman 

(14),    Gruber  and  Krupke  (IS),    and  Gruber  et  al.    (16). 

169 
Tm  appeared  to  be  an  isotope  particularly  suited  for  studies 


(13)  A  preliminary  report  of  part  of  this  work  appeared  elsewhere 
(7). 

(14)  E.  Y.    Wong  and  I.    Richman,    J.Chem.    Phys.  J^    1182(1961) 

(15)  J.B.    Gruber  and  W.F.    Krupke,    to  be  published 

(16)  J.B.    Gruber,    W.F.    Krupke  and  J.  M.    Poindexter,    to  be 
published 


-6- 
of  electronic  shielding,    for  the  following  major   reasons: 

(1)  The   low  energy  of  the   8.4   keV  transition  used  results   in  a  high 
Debye-Waller  factor  (recoil-free  fraction)  even  at  very  high  tempera- 
tures,   thus  permitting  a  measurement  of  the  quadrupole  interaction 
within  an  unusually  wide  temperature  range. 

(2)  The   separation  of  the  excited  levels   belonging  to  the  ground 
multiplet  of  thulium  (L  =  6;  S  =    I)  is   rather  large,    with  the  first 
excited  level  (    H.)  some   5600  cm        above  the  ground  term  (    H/). 
Thus  the  existence  of  the  higher   levels  of  the  ground  multiplet  is  of 
minor  concern  for  the  interpretation  of  our  data  in  thulium,    in  con- 
trast to  the  situation  prevailing  in  the  case  of  some  other  rare  earth 
ions. 

(3)  The  spin  of  the  nuclear  ground  state  (I  =  1  /2)  and  of  the  8.4  keV 
excited  state  (1  =  3/2)  is  rather  low  ,  resulting  in  a  small  number  of 
quadrupole  hyperfine  components  of  the  gamma  lines  which  are  easily 

resolvable. 
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(4)  The  nuclear  collective  model  applies  well  to  Tm  thus  per- 
mitting a  rather  reliable  semi -theoretical  estimate  of  the  nuclear 
quadrupole  moment  of  the  8.4  keV    state. 

(5)  The  relative  abundance  of  Tm  is    100%. 


-7- 
II.      CRYSTAL   ELECTRIC    FIELD   (CEF)    INTERACTIONS 

A  rare  earth  ion  interacts   in  a   salt  with  the  CEF  produced  by  all 
the  ions  which  surround  its  position  in  the  lattice.     The  dominant 
effect  is  the  interaction  of  the  CEF  with  the  electrons  in  the  partially 
filled  4f-shell.      This  interaction  is  weak  compared  to  the  spin-orbit 
interaction,    in  contrast  with  the  situation  prevailing  in  the  case  of 
iron-transition  elements.     As  a  result,    the  total  angular  momentum 
J    remains  a  good  quantum  number  for   rare  earth  ions  bound  in 
crystals.      The  effect  of  the  CEF  then  essentially  is  a  partial  or  com- 
plete removal  of  the  ZJ  +    1  fold  spatial  degeneracy  of  the  orientation 
of   J    which  exists  in  a  free  ion.      The  actual  number  of  electronic  CEF 
levels  depends  on  the  symmetry  of  the  field,    while  the  level  spacing 
depends  on  the  strength  of  the  interactions  between  the  CEF  and  the 
4f  electrons.     The  situation  is  illustrated  in  Fig.     la. 

The  potential  energy  describing  the  interaction  between  the  CEF 
and  a  negative  charge  at  position  (r,   d  ,  v?  )  within  the  ion  centered  at 
the  origin  may  be  represented  in  good  approximation  by  the  following 
expansion,    not  including  shielding  from  closed  shells: 

+  n 

-eV(r,^,  y  )  =     y        T]       Am   rn    <ZP    (^,   «/)       ,  (1) 

<t_ i        ■     *  n  n 

n       m  =  — n 

if  one  assumes  that  there  is  no  overlap  between  the  charge  distri- 
butions of  different  ions.     In  Eq.    (1)  the  A      represent  lattice  sums 
over  point  charges  and  effective  multipole  moments  in  the  surround  - 
ing  ions.      The  relevant  functions    0       ,    which  are  linear  combinations 
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Fig.    la:         Schematic  of  the  atomic  level  splitting  of  a  rare  earth 

ion     in  the  CEF.      For  a  nuclear  spin  1=3/2  the  nuclear 
quadrupole  interaction  splits   each  CEF  level  into  a 
doublet,    which  is  the  case  illustrated.      Typical  overall 
CEF  splittings   are  of  the  order  of   10    "  eV,    while  typical 
quadrupole  hyperfine   splittings  are  of  the  order  of   10 
eV. 

Fig.    lb:         Schematic  of  the  nuclear  quadrupole  splitting  of  the 
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8.  4  keV  transition  in  Tm         .      The  temperature  dependent 

level  splitting  (AE)T ,    which  is  typically  of  the  order  of 

10        eV,    is  the  average  of  the  hyperfine  splittings  of 

Fig.    la,    weighted  according  to  their  Boltzman  factors. 
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of  spherical  harmonics   Y         and    Y        ,    are  defined  as  follows  (17) 

n  n 


2n 


.  .2n)P,     (cos^) 
en 


A±m       ,  ,,m  .      ,.   (n  — m)   '.        ,-.m   .  „q,. 

0  =     2  rn'.f- T — Vt"      P         cos  XT)     • 

n  (n+m).  n 


[cos  m^ 
'sin  m*/> 


where  P      and   P       are  Legendre  polynomials  and  associated  Legendre 
functions,    respectively.     In  particular,    we  obtain  for    n  =  2: 


0°z  =  3  cos2  #-    1 


(2a) 


0CZ  =  sin2    ^cosZ/7 


(2b) 


0"Z=   sinz  1?  sin  2  Z7 


(2c) 


Specifically,    the  Hamiltonian  describing  the  interaction  between 
the  CEF  and  the  electrons  in  the  partially  filled  4f  shell  of  rare  earth 
ions,    including  the  effect  of  shielding  via  the  closed  electron  shells 
of  the  central  ion  is   given  by: 


H 


(4f) 

CEF 


k 


A 


m 


n  m 


k  n       k 


«£<t>k.*k>     (3: 


(17)  The  normalization  of  the  functions  0       {rv,  't  )  is  arbitrary; 

the  choice  adopted  here  is  the  one  most  commonly  used 
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The  terms  proportional  to    r       describe  the  potential  energy  due  to 
the   direct  interaction  of  the   CEF  with  the  k-th  electron  in  the  4f-shell 

while  the  terms  proportional  to   S    (r.  )  describe  the  additional  poten- 

r      r  n      k  ' 

tial  energy  arising  from  a  deformation  of  the  closed  electron  shells. 
The  interaction  described  by  the  Hamiltonian  in  Eq.(3)  splits  the 
electronic  ground  state  of  the   free  ion,    characterized  by  total  angu- 
lar momentum    J    ,    into  a  number  of  CEF  levels.      We  shall  assume 
in  calculating  these  CEF  levels  that  the  angular  and  radial  parts   of 
the  free  ion  wave  functions  can  be  factorized  and  that  higher  terms 
with  different    J    values  can  be  neglected.      Under  these  circumstances 
we  are  dealing  with  a  manifold  of  states   belonging  to  the  same    J    and 
it  is  then  convenient  to  replace  the  angular  operators  occurring  in 
the  Hamiltonian,    Eq.(3),    by  equivalent  operators  (18).     The  relevant 
matrix  elements  then  are  of  the  form 

H  =  y     Am    /rn>„     <J||d   ||  J><J,mT|om(J    ,J    ,J    )|j,m'>    (4) 

mTmTI       Z-»         n      \     'E      x    "     n"     '     >  Jnx       y      z1  J 

J      J '       n  m 

where  ^r^   =  ( 1  ^    <rn)4f  (5) 


(18)  K.W.H.    Stevens,    Proc.    Phys.    Soc.    A6_5,    209(195Z); 

R.J.  Elliott  and  K.  W.H.  Stevens,  Proc.  Roy.  Soc.  A2 1 8  , 
553  (1953);  J.  P.  Elliott,  B.  R.  JuddandW.A.  Runciman  , 
ibid.    A240,    509(1957);     R.    Orbach,    ibid.    A264  ,    458(1961) 
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n 


<« 


4f 


Sn(r) 


U4f>/^n>4r  <6> 


X    '4f  \    4f  I      nl        4f/ 


(7) 


and 


t> 


a,    3-    7        for        n  =  2,    4,    6 


In  these  expressions  U  .,   is  the  radial  part  of  the  electronic  wave 


4f 


.m, 


functions  for  the  4f-shell.      The  functions  O      (J    ,  J    ,  J    )  are  opera- 

n       x       y       z  r 

tor  equivalents;     those   relevant  for  this  work  are   listed  in  Table  I. 
The  expressions       \J       ru     ||J  /    are   reduced  matrix  elements  (19), 
which  for    the  more  general  case  of  intermediate  coupling  are  avail- 
able for  Tm         in  the  literature  (  14)-(  17),  (20). 


m 


It  is  in  principle  possible  to  calculate  the  parameters  A         and 


n 


\r      A^,    but  difficult  in  practice.      Difficulties  are  in  the  evaluation 


m 


of  the  "lattice  sums"  A       because  of  a  lack  of  sufficient  knowledge  of 

n  ° 

the  ionic  position  coordinates  and  their  temperature  dependence  as 
well  as  of  the  values  of  moments  in  the  surrounding  ions  (21).      The 
evaluation  of  the  radial  integrals     \r    /    ,    which  are  the  expectation 
values  of     r       for  the  4f  shell  modified  by  contributions  from  closed 
shells  to  the  electric  multipole  fields  at  the  4f  electron  positions,    is 
hampered  by  the  lack  of  knowledge  of  sufficiently  accurate  atomic 


(19)  B.R.    Judd,    Proc.    Roy.    Soc.    A24 1 ,    414(1957) 

(20)  J.  B.    Gruber  and  J.G.    Conway,    J.    Chem.    Phys.    32   ,     1531 
(I960) 


(21)  M.T.    Hutchings  and  D.K.    Ray,    Proc.    Phys.    Soc.    81,663 

(1963) 
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wave  functions  for  bound  rare  earth  ions.      For  these   reasons  it  is 
therefore  preferable  to  introduce  the  "CEF  parameters" 


~m  Am     /   ri\  ,  Q, 

Cn       -"     An      V   /E  (8) 


to  be  determined  by  experiment.      The  point  symmetry  of  the  central 
ion  drastically  reduces  the  number  of  CEF  parameters  (22).      In  the 
case  of  rare  earth  ions  only  the  terms  with    n  =  2,  4,    6    need  to  be 
considered,    with  the  effects  of    n  =   1,  3,  5    being  negligible  in  most 
cases  (22). 

The  wave  functions      r\.    of  the    V*-th  CEF  level  will  be  taken  as  a 
linear  combination  of  eigenvectors  of  the  total  angular  momentum    J  , 


(mA,mj 


%-    U4f<"    E-^V  <9> 


y  mj 


(«.,) 


The  expansion  coefficients     a    .    '         and  the  energy  eigenvalues   E  j 
follow  from  the  diagonalization  of  the  interaction  matrix  H 
Eq.   (4). 


m,    m' 

■J    ,  %J 


(22)  A  compilation  of  the   relevant  values     n     and     m     for  various 

crystal  symmetries  was  given  by  J.  L..    Prather,    NBS 
Monograph    19  (1961) 
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III.     THE    NUCLEAR    QUADRUPOLE   INTERACTION 
Each  of  the  CEF  levels  may  produce  a  magnetic  field  and  an 
electric  field  gradient  at  the  nuclear  site;    this  results  in  hyperfine 
splittings  of  the  electronic  levels.     A  rare  earth  nucleus  thus  ex- 
periences at  a  certain  time  a  magnetic  field  and  an  electric  field 
gradient  which  depends  on  the  electronic  state  that  is  actually  popu- 
lated at  this  time.    The  situation  substantially  simplifies  at  elevated 
temperatures  where  the  spin-lattice  relaxation  phenomenon  produces 
rapid  transition  between  the  different    CEF  levels.      The  nucleus 
under  these  circumstances  experiences  a  magnetic  field  and  an  elec- 
tric field  gradient  which  in  essence  result  from  averaging  these 
fields  over  all  electronic  states  weighted  according  to  the  population 
numbers.      This  averaging  process,    which  essentially  constitutes  a 
time  averaging  process,    holds  only  if  the  significant  electron  relax- 
ation times  are   short  compared  to  all  other   relevant  times   such  as 
the  nuclear  lifetimes  and  the  nuclear  precession  times,    a  situation 
prevailing  at  temperatures  above  a  few  degrees  Kelvin.     In  particu- 
lar,   the  magnetic  hyperfine  interaction  cancels  in  the  absence  of  an 
external  magnetic  field  and  all  one  is  left    with  is  the  quadrupole 
hyperfine  interaction  (23),  (24).      An  example  of  this   situation  is 
illustrated  in  Fig.    lb  for  an  assembly  of  nuclei.      The  quadrupole 


(23)  R.  L.    Cohen,    U.    Hauser  and  R.  L.    Mdssbauer,    Proc. 
MHssbauer  Conf.    2nd,     (John  Wiley  and  Sons ,    N.  Y.  ,  1962) 
p.  172 

(24)  R.  L.    Mdssbauer,    Rev.    Mod.    Phys.    36,    362(1964) 
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interaction  is   strongly  temperature  dependent  since  the  overall  CEF 
splitting  within  the   lowest  electronic  term  is   only  of  the   order   of  a 
few  hundred  cm 

The  electric  field  gradient  which  interacts  with  the  nuclear 
quadrupole  moment  of  a  rare  earth  nucleus  bound  in  an  ionic  crystal 
has  four  significant  sources: 

1)  One  contribution  is  the  direct  field  gradient  produced  at  the 
nuclear  site  by  all  of  the  ions   surrounding  the  host  ion  which 
c  ontains  the  nucleus  in  question.      This   contribution  in  the 
case  of  rare  earth  ions  is  usually  negligible  in  comparison 
with  the  contributions  from  other  sources,    particularly  at 
low  temperatures. 

2)  Another  contribution  results  from  the  electric  field  gradient 
produced  at  the  nuclear  site  by  the  electrons  in  the  partially 
filled  4f-shell  of  the  host  ion.      This   field  gradient  results 
from  the  interaction  of  the  4f-electrons  with  the  CEF  pro- 
duced by  the  surrounding  ions.      This  interaction  effectively 
induces  electric  multipole  moments  (multipole  polarization) 
in  the  4f-shell;     the  quadrupole  part  of  this  polarization  con- 
tributes to  the  electric  field  gradient  experienced  by  the 
nucleus . 

3)  A  distortion  is  usually  also  induced  by  the  CEF  in  the  closed 
electron  shells,    yielding  another  contribution  to  the  total 
field  gradient  experienced  by  the  nucleus.      This  contribution 
is  proportional  to  source  (1),    with  proportionality  factor   -y 
The  absolute  value  of  the  proportionality  factor  is  in  the 
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case  of  the  rare  earths  usually  large  in  comparison  with 
unity,    thereby  leading  to  such  an  enhanced  field  gradient 
(antishielding  effect)  that  it  often  becomes  comparable  with 
the  one   resulting  from   source  2).      This  is  the   "lattice" 
Sternheimer  effect  (25)-(28). 
4)  Another  field  gradient  contribution  due  to  an  induced  quadru- 

pole  moment  in  the  closed  electron  shells   results  from  the 
interaction  of  the  closed  electron  shells  with  the  electrons  in 
the  partially  filled  4f-shell.      This   relatively  small  contri- 
bution,   which  is  proportional  to  source  2),    with  proportionali- 
ty factor   —  R/~v.    is  the  "atomic"  Sternheimer  effect  (5),  (28). 
Collecting  the  different  contributions,    we  obtain  for  any  compo- 
nent   eq  .      of  the  electric  field  gradient  tensor 


ecl,i    ^d-Yj^0    +     ^-R0)e^i]       i.j  =  1.2.3  (10) 


where  y        and    R~    are  the  lattice  and  atomic  Sternheimer  factors, 
'co  U 

respectively,    as  introduced  above. 


(25)  E.G.    Wikner  and  G.    Burns,    Phys.    Letters_2_,    225(1962) 

(26)  D.K.    Ray,    Proc.    Phys.    Soc._82,    47(1963) 

(27)  R.M.    Sternheimer,    Phys.    Rev.    132,    1637(1963) 

(28)  A.J.    Freeman  and  R.  E.    Watson,    Phys.    Rev.    13_2  ,    706 
(1963) 
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In  the  principle  axes  system  of  the  electric  field  gradient  tensor 

(V  ) 
the  nuclear  quadrupole  interaction  Hamiltonian  H*  associated  with 

the  (V)-th  CEF  level  of  the  ion  is  given  by 

(11) 


where    I  ,    I.    ,  J[       are  the  usual  nuclear  spin  operators  and  Q  is  the 
nuclear  quadrupole  moment.      The  quantities     q. .  and\y|^..       VO 

determine  the  direct  contributions  to  the  electric  field  gradient  at 
the  nuclear  site  produced  by  the  surrounding  ions  in  the  lattice  and 
by  the  4f-electrons  of  the  host  ion,    respectively.     They  are  defined 
by 

eq(£at)=   [^Vfr^^l/h.       B  x.  ]  r  =  Q     ,  (12) 

where     V(r,/2r,  V)  is  defined  in  Eq.  (1),    and 

(13) 
....     .  .       4f-electrons  ,        .  x        / 


where  the  wave  function    |  \} \    of  the     \l-th  CEF  level  is  of  the  form 
given  by  Eq.  (9). 

Explicitly  we  obtain  for  the  lattice  contribution  from    Eqs.(l), 

(2), (12) 


2       (Lat)  A  A  0  z,    (Eat)         (Eat).  ..2  .... 

ez    qv         '=-4A  ;      eMq  -q  ')=-4A  14 

nzz  2  xx  yy        '  2 
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Expressing  the  contributions  froin  the  4f  -electrons  within  a  manifold 
of  states  of  constant     J     in  terms   of  operator  equivalents,    we  obtain 
from  Eq.  ( 13) 

<v|alt>>    =-<JH°«J>   <*'\  <"|3£-j,,|»'>  <15*> 

(15b) 

where     \r      y  f  is  defined  by  Eq.(7). 

Usually  one  observes  only  an  average  field  gradient  from  the  4f- 
electrons,    which  is  a  field  gradient  from  the  individual  CEF  levels 
weighted  according  to  their  Boltzman  factors,    as  discussed  above. 
If  we  consider  only  those  electronic  states  which  belong  to  the  lowest 
manifold  spanned  by  the  state  vector  J^  ,    then  the  average  direct 
contribution  from  the  4f-electrons  to  the  electric  field  gradient  acting 
on  the  nucleus  at  temperature     T     is  given  by  (29) 


(4f)  „.   , 


2J  +  1 
£      <^U!ffV>     •     exp(-E»//kT) 

<qii>T=    — (i6) 

2J  4-  1 
£        exp(-Ev/kT) 

V=   1 


The  diagonal  component  of  the  averaged  total  electric  field  gradient 
tensor  is  according  to  Eqs.(10)  and  (16)  given  by 


(29)  A  more  general  description  including  effects  of  higher      J 

states  is  given  in  the  Appendix  I 


•  21 


/  v         -    /i  \     (Lat) 


(if) 


OO  11 


+  (1-V<s!i1"> 


(W) 


where  we  have  neglected  any  temperature  dependence  of  the  lattice 

contribution  q.. 
ni 

The  total  Hamiltonian  describing  the  average  quadrupole  inter- 
action at  temperature  T  may  now  according  to  Eqs.  (  1  1)  and  (  17)  be 
written  as 

(18) 
e2Q 


HQ(T)  - 


41(21-1) 


<0  t  <3r,-l2» +  <2xx-ayy>T  K£i  +i2-» 


We  shall  now  apply  the  preceding  formalism  to  the  particular 

169 
case  of  Tm         .      The  twofold  degeneracy  of  the  nuclear  ground  state 

1  69 
of  Tm  (1=   l/2)is  not  removed  by  the  Hamiltonian,    Eq.  (  1  8) ;     the 

8.4  keV  excited  state  (I  =  3/2),    on  the  other  hand,    is  split  by  the 
nuclear  quadrupole  interaction  into  two  states.      Their  energy  separa- 
tion/aE")  T    which  follows  from  the  diagonalization  of  the  Hamiltonian 
Hq(T)  is  given  by 


/aE>„  =  (e2Q/2)      <q       >T     +    |    <q       -q      \ 


yy/T 


(19) 


This  expression  may  be  written  in  more  detail,    by  using  Eqs.  (5),    (8), 
(14), (15), (17): 
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.0 


^Me»Q/2){    [<J[.iJ>    <r-3>Q<3^-J^T+-^.      {l.yJ 


<<■   . 


f 


4C 


2 


i     1 

2 


^|  OII«||J><r-3>Q<j2++j2->T+^T7    (1-^co>J      } 

E  (20) 

where   <^3J2   -  J2/_,     and      O2     +   J2  Y_     are  thermal  averages  defined 
as  those  given  by  Eq.(l6),    while  the  parameter    (r       /-.     is  defined  by 

<r-3>Q  =  (l-RQ)      <r"3>4f  (21) 

It  is  just  this   splitting   (AE)T    that  is  measured  as  a  separation  of 
gamma  lines  in  recoilless   resonance  absorption  experiments. 

Several  additional  hyperfine  interaction  mechanisms  which 
contribute  to  the  net  nuclear  quadrupole   coupling  of  a   rare  earth  ion 
have  been  neglected  in  our  calculations.      These  additional  contri- 
butions arise  in  second-order  perturbation  theory  with  the  principal 
effects   conning  from  the  magnetic  hyperfine  interaction  itself  (30) 
(the  so-called  pseudo-quadrupole  coupling)  and  from  the  quadrupole 
interaction  with  states   of  higher     J     admixed  into  the  ground  state 
multiplet  by  the  CEF.      We  have  made  calculations  of  these  contri- 
butions for  the  compounds  covered  in  this  paper  and  they  amount  to 
less  than   1%    of  the  total  quadrupole  interaction  energy. 

In  order  to  compare  experimental  results  with  theory  within  the 
framework  of  the  CEF  model  it  is   convenient  to  replace  in  the  theo- 


(30)  R.J.    Elliott,    Proc.    Phys.    Soc.    B70  ,    119(1957) 
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retical  expression  for  the  quadrupole  splitting    <(AE^>T  all  quantities 
involving  electronic   radial  integrals  (the  theoretical  determinations 
of  which  is  presently  somewhat  uncertain)  as  well  as  the  nuclear 
quadrupole  moment  by  experimentally  observable  parameters.     For 
this  purpose  we  introduce  the  dimensionless  parameters 


pl   =e 


2Q     <r"3>0    <j||clJ>    /C? 


'Q 


(22a) 


=   Q(i-y(X))/<r2> 


(22b) 


Expressed  in  terms  of  these  parameters  the  quadrupole   splitting 

169 
in  Tm  reduces  to 


<AE>T=    i{     [c°p,     <^Z-LZ\     +4C°    P, 


(23) 


4 


(3/2>q    p,    <  J+2    +   J_2>T 


+  4CI  Pz 


The  temperature  averages    <  3jJ  -  J.2^     and    <J_^    +    J^2_  >T    within 

the  framework  of  the  CEF  model  depend  only  on    the  experimentally 

m 
observable  CEF  parameters     C 


n 
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IV       EXPERIMENTAL   TECHNIQUE 
The  nuclear  quadrupole  interaction  was   measured  by  using  the 

technique  of  the  recoilless  nuclear  resonance  absorption  of  gamma 

169 
radiation  (2).      The  partial  decay  scheme  of    Er  is   shown  in  Fig.   2. 

Measurements  of  the  gamma   resonance  absorption  were  performed 

as  a  function  of  the   relative  velocity  between  sources  and  absorbers. 

The  measurements  involved  sources  of  erbium  trifluoride  (ErF,)  and 

erbium  oxide  (Er.O,)  and  absorbers  of  thulium  ethyl  sulfate 

(Tm  (C-H-SO.)-    •    9H..O  ,    abbreviated  to  TmES)and  thulium  oxide 
L     d       4   3  L 

(Tmz03). 

Anhydrous   ErF,.   provides  an  excellent  source  for  experiments 

169 
utilizing  the   8.4  keV  line  of  Tm         .      The  crystal  structure  of  the 

heavy  rare-earth  trifluorides  has  been  investigated  by  Zalkin  and 
Templeton  (3  1).      At  temperatures  below  about  900  -  1000°  C  the 
stable  phase  is   orthorhombic  ,    space  group  D        -Pnma,    having  four 
formula  units  per  unit  cell.      The  rare  earth  ions  are  c  rystallogra  - 
phically  equivalent,    having  the  point  symmetry    m.      Thus,    although 
the  electric  field  gradient  tensor     eq  .      is   not  axially  symmetric,    all 
erbium  (or  thulium)  nuclei  experience  the  same    eq...       Therefore, 
the  quadrupole  splitting  of  the   recoilless  absorption  line   given  by 
Eq.(23),    may  be  expected  to  pass  through  zero  or  at  least  through  a 
minimum  at  a  specific  temperature  (550    K  in  this  case).      The  advan- 
tages of  a  single-line   source  are  thereby  obtained.      The  line  width 
obtained  this  way  with  sources  of  ErF,   is   less  than  with  sources  of 


(31)  A.    Zalkin  and  D.  H.    Templeton,    J.    Am.    Chem.    Soc.    75_, 

2453  (1953) 
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Fig.    2:  Partial  decay  scheme  of    Er 
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Er.,0     where  the  presence  of  non-equivalent  erbium  sites  complicates 
the   situation  (23),  (24).      At  the   same  time,    with  reasonable  precau- 
tions,   the  ErF     can  be  maintained  at  the  critical  temperature  for 
periods  of  several  weeks  without  decomposition  or  reaction.      Tins 
chemical  stability  does  not  exist  with  most  other  erbium  salts  in 
which  the  erbium  ions  are  also  crystallographically  equivalent  (e.  g.  , 
the  sulfate,    nitrate,    chloride). 

Anhydrous   ErF,  was  prepared  from  erbium  metal  or  erbium 
oxide  by  a  "wet"  process.      The  metal  or  oxide  was  first  dissolved  in 
a  small  quantity  of  nitric  or  hydrochloric  acid  in    polyethylene  centri- 
fuge tube.     A  few  ml  of  aqueous  hydrofluoric  acid  were  then  added 
and  the  mixture  heated  at  approximately   100    C  in  a  water  bath  for  30 
minutes.      The  somewhat  gelatinous   ErF,   precipitate  was  then  centri- 
fuged  down,    the  excess  solution  decanted  off,    the  precipitate  washed 
with  distilled  water,    centrifuged  three  to  five  times  and  dried  in  air 
at  roughly   100     C.      Air  drying  yields  a  hydrated  ErF,   of  unknown 
composition.      To  remove  the  water  of  hydration,    the  dry  contents   of 
the  centrifuge  tube  bottom  were  transferred  to  a  small  tantalum  boat 
and  annealed  in  an  evacuated  fused  quartz  tube.      Experience  showed 
that  the  hydrated  ErF     could  be  converted  directly  into  a  mixture  of 
the  several  forms   of  oxyfluoride  (32)  if  the  annealing  temperature 
was   raised  too  rapidly.      The  procedure  finally  adopted  was  to  hold 
the  hydrate    at  room  temperature  at  about   10        torr  for  at  least   12 
hours  in  order  to  pump  off  most  of  the  water.      The  temperature  was 


(32)         W.  H.    Zachariasen,    Acta  Cryst.    4,    231  (1951) 
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then  raised  slowly  (in  6  hours)  to   150     C,    thus    removing  virtually  all 
of  the  water.      Finally,    the  temperature  was   raised  to  850     C   in  an- 
other 6  hours  and  then  reduced  back  to  room  temperature  within  I 
hours.      This  procedure  yielded  consistently  good  clean  x-ray  powder 
patterns  of  the  orthorhombic  phase  without  a  trace  of  the  hexagonal 
phase  appearing  (31).      ErF^  prepared  in  this  manner  appears   to   re- 
main stable  at  room  temperature  over  an  indefinite  period  of  time. 
At  elevated  temperatures  care  must  be  ex^ercised  to  avoid  reaction 

with  oxygen  or  water  vapor.      Sources  of  ErF..   were  prepared  in  the 

1  ft  H 
above  manner  from  Er,0^  (usually  enriched  in  Er         )  or  from  erbium 

metal  after  irradiation  in  the  Materials  Testing  Reactor,  Arco,  Idaho. 
Alternatively,  the  ErF,  was  prepared  first  and  then  irradiated.  Iden- 
tical spectra  were  obtained  by  the  two  methods. 

Absorbers  of  TmF,  were  used  in  order  to  experimentally  determine 
the  critical  temperature  at  which  the  narrowest  possible  emission 
line  is  obtained  with  sources  of  ErF...      Figure   3  shows  the  tempera- 
ture dependence  of  the  quadrupole  splitting  in  TmF.,.      The   source    was 
mounted  in  a  small  evacuated  oven  shown  in  Fig.    4.      The  absorber 
was  maintained  in  a  helium  atmosphere  within  an  oven  equipped  with 
beryllium  windows. 

It  is  interesting  to  note  that  the  same  minimum  line  width 
(1.8  cm/sec)  was  obtained  in  both  the  trifluoride-trifluoride  and  tri- 
fluoride-ethylsulfate  experiments.      This   strongly  suggests  that  the 
quadrupole  splitting  of  the  trifluoride  line  does  indeed  pass  very  near 
to  zero  at  550°  K  (24).      This  minimum  observed  line  width  of 
1.  8  cm/sec  may  be  compared  with  the  theoretically  predicted  line 
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Fig.    3:  Temperature  dependence  of  the  quadrupolc  splitting 

169 
of  the   8.4  keV   gamma  line  of  Tm  using  an  ErF, 

source  and  a  TmF.  absorber.     Source  and  absorber 

were  maintained  at  the  same  temperature. 


29- 


-24.0 


200  400 

TEMPERATURE    (°K) 


FIG.    3 


-30- 
Fig.    4:  Details  of  moveable  source    oven.      The  entire  oven 

(weight   0.  2Kg)  was   moved  relative  to  the   absorber 
by  the  cam  drive.      The  main  body  of  the  oven  was 
made  of  stainless   steel.      A  similar  oven  made  of 
aluminum  was  used  with  the  transducer  drive.      The 
heating  element  was  fabricated  from  nichrome  strips, 
1/16  x  0.  005  in.      For  a  source    temperature  of  550     K 
the  power  input  was  30  W. 
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width.     In  a  transmission- vs -velocity  measurement  since  an  emission 
line  of  width     T    is  moved  over  an  absorption  line  of  width    T   ,    one  ex- 
pects a  minimum  line  width  of  2T   .     In  the  case  at  hand  this  is  (in 
velocity  units) 


2Tc/E      =     2'fic/'CE      =     0.74  cm/sec 

y  y 


-9 
based  on  the  lifetime  (33),  7.    ~     6. 28*    10        sec,    and  energy, 

169 
E      _    8.42  keV,    of  the  first  excited  state  of  Tm         .      However  for  an 

y 

absorber  of  finite  thickness  this  width  increases  by  a  correction 
factor  which  is    1.47  in  our  case  for  an  absorber  of  5  mg/cm2  of 
thulium  and  a  total  conversion  coefficient  of  325  from  Kankeleit  et  al. 
(34).      Thus  the  predicted  line  width  is    1.09  cm/sec  which  must  be 
compared  with  the  minimum  observed  line  width  of  1 .  8  cm/sec.      The 
observed  line  width  is   1.  6  times  broader  than  expected.      The  origin 
of  this  line  broadening  is  uncertain. 

Absorbers   of  TmES  were  prepared  by  crushing  single  crystals. 
Absorbers  of  Tm?0_  and  sources  of  (enriched)  Er    O,  were  prepared 
from  commercially  available  material.     Absorbers  of  all  materials 
to  be  used  below  room  temperature  were  prepared  by  mixing  the 
powdered  samples  with  a  soft  wax  and  pressing  the  mixture  into  thin 
disks  between  mylar  films.     Absorbers  and  sources  of  all  materials 


(33)  R.  E.    McAdams,    G.  W.    Eakins,    E.N.    Hatch,    Phys.    Letters 
6,    219  (1963) 

(34)  E.    Kankeleit,    F.    Boehm,    R.    Hager,    Phys.    Rev.    in  press 
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to  be  used  above   room  temperature  were  prepared  by  settling  the 
powdered  samples   from  a  slurry  of  dry  acetone  onto    1 /Z  mm  thick 
beryllium  windows. 

The   relative  velocities   required  for  Doppler -shifting  the  gamma 
lines  were  produced  by  using  both  cam  drives  (35)  (providing  constant 
velocities)  and  transducer  drives  (  36)  (providing  constant   acceleration). 
A  block  diagram  showing  the  experimental  apparatus  for  use  with  the 
cam  drive    is   shown  in  Fig.    5.      The  experimental  arrangement  used 
with  the  transducer  drive    is  given  in  Ref.(36).      Proportional  counters 
filled  with  one  atmosphere  of  a  mixture  of  90%  argon  and   10%  methane 
(by  volume)  and  equipped  with   1 /Z  mm  thick  beryllium  windows  were 
used  as   detectors,    see  Fig.    7. 

A  cryostat  specifically  designed  for  recoilless   resonance  absorp- 
tion experiments  with  low  energy  gamma  radiation  was  used  for  the 
measurements  (37).      The   sample  temperatures   in  the  range  from    10    K 
up  to   300     K  were  attained  by  either  controlled  heating  of  the  cooled 
sample  holder,    by  pumping  on  liquified  gases,    or  by  using  exchange 
gas  cooling.      The   sample  disks  were  clamped  between  thin  beryllium 
disks   soldered  to  the  cryostat  sample  holder  in  order  to  ensure  good 
temperature  uniformity  and  stability.      Temperature  measurements 
were  made  using  carbon  resistors  and  thermocouples.      The  oven  used 

for  heating  sources  to  ZO00     K  is   shown  in  Fig.    8. 

1  69 
Some  typical  Mbssbauer  spectra  for  Tm  are   shown  in  Fig.    9. 


(35)  R.  L.    Mtissbauer,    Proc.    MOssbauer  Conf.    Znd,  (John  Wiley 
and  Sons,    New  York,     1962)  ,    p.  38 

(36)  E.    Kankeleit,    Rev.    Sci.    Instr.    3_5,    194(1964) 

(37)  F.    T.    Snively,    to  be  published 
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Fig.    5:  Block  diagram  of  experimental  apparatus  for  use  with 

the  cam  drive.      A  detailed  schematic   of  the  pro- 
grammer  is  presented  in  Fig.    6.      All  other  electronic 
equipment  shown  is   commercially  available.      The  cam 
drive  has  been  described  elsewhere  (35).      The    recoil- 
less   resonance  absorption  measurements  were  per- 
formed by  first  moving  the  source  at  a  normalizing 
speed  (12  cm/sec),    then  at  a  measuring   speed  (v),    and 
finally  at  the  normalizing  speed.      The  three   runs  were 
of  about   5  min  each.      The  counting  rates,    C,    during 
the  two  normalizing  runs  were  averaged  and  combined 
with  the   results  of  the  measuring  speed  to  yield  the 
amount  of  absorption,    A(v). 

A(±v)  -  [C(±12)AVE-C(±v)]/C(±12)AVE 
This   sequence  of  events  was   repeated  until  the  counting 
statistics  were  satisfactory.      The  information  at  the  end 
of  each   run  was  printed  out  on  a  typewriter   for  moni- 
toring purposes  and  punched  out  on  paper  tape  for 
processing  by  an  electronic  computer.      The  entire 
process  was  automatic  except  for  changing  from  one 
measuring  speed  to  another. 
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-36- 
Fig.    6:  Schematic  of  programmer  for  use  with  cam  drive. 

This   programmer   controlled  the   changes   from   the 
normalizing  speed  to  the  measuring  speeds  and 
initiated  the  counting  and  print-out  cycles  of  the 
scaler,    see  Fig.    5. 
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FIG.    6. 
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Fig.7:  Proportional     counter.      The  anode   of  the  counter  is 

a  3  mil  stainless   steel  wire  kept  under  light  tension 
by  a  spring  in  the  end  fitting.      The  beryllium  window 
is   0.  5  mm  thick.      This   counter  was   operated  at  about 
2500  V  with    1   atm  of  90%  argon  and   10%  methane,    by 
volume.      Under  these  conditions  the   resolution  at 
8.42  keV  was  about   17%. 
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Fig.    8:  High  temperature,    resistance  heating,    source  oven. 

The  heating  coil  was  fabricated  from  two  pieces  of 
56  mil  tungsten  wire   15  in  long.      These  two  pieces 
were  connected  in  parallel  and  the  source  container 
was   suspended  from  the  midpoint  of  each  piece.      For 

a  source  temperature  of   1970     K  the  power  input  was 

_4 
4.  8  kW.      A  vacuum  of   10        torr  was  maintained  in 

the  oven.     The  thermocouple  made  of  tungsten  vs. 

tungsten  -  26%  rhenium  is  reliable  to  3100    K. 
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169 
Fig.    9-'  Quadrupole  splitting  of  the   8.  4  keV  level  of  Tm 

in  an  absorber  of  thulium  ethyl  sulfate  (5  mg/cm2 

169 
of  thulium).     A  "single  line"  source  of  Er  in  ErF, 

was  used  at  the  critical  temperature  T  =  550     K 

throughout  curves  a-d.      The  spectra    a,    b    and    c,    d 

were  obtained  by  using  a  cam  drive  and  a  transducer 

drive,    respectively. 
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V.      EXPERIMENTAL    RESULTS   AND   ANALYSIS 

The  nuclear  quadrupole  interactions  measured  as  a  function  of 
temperature  in  the  compounds  TmES  and  Tm,0,  are  given  in  Figs. 
10  and   11,    respectively.      Details  of  the  figures  are  explained  below. 
The  reduction  of  our  experimental  results  is  carried  out  in  two  sub- 
stantially different  ways: 

Method   1:     We  combine  our  nuclear  quadrupole  splittings  obtained 
from  gamma  resonance  absorption  measurements  with  optically  de- 
termined CEF  levels  and  obtain  two  quantities 

PlC°   [(C^)2    +     ]-    (C^V      and        PZ[(C°)2+      i-(C^)2]1      which  de- 
pend  directly  on  the  electronic  shielding  factors,    compare  Eqs.  (5), 
(21),    (ZZ).      This  method  emphasizes  the  low  temperature  data,    which 
have  the  smallest  relative  errors. 

Method  Z:     The  same  two  quantities  may  be  obtained  without  the 
necessity  of  referring  to  any  optical  determination  of  CEF  levels, 

merely  by  using  gamma  resonance  measurements  obtained  at  elevated 

h 

temperatures.      This  method  is  useful  in  those  cases  were  measure- 
ments can  be  performed  at  temperatures  which  are  large  compared 
to  the  overall  CEF  level  splitting,    but  small  compared  to  the  spin- 
orbit  splitting.     This  is  the  case  in  both  TmES  and  Tm?0~. 

1.      Thulium  Ethyl  Sulfate     (TmES) 
All  rare  earth  lattice  sites  in  TmES  are  occupied  by  Tm         ions 
with  point  group  symmetry  C       .      By  choosing  the  proper  coordinate 
system  (ZZ)  the  relevant  CEF  parameters  as  defined  by  Eq.  (8)  are 
limited  to  C?,    C.,    C,     and  C,   for  this  symmetry.     This  leads  to  an 
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Fig.    10:  Temperature  dependence  of  the  nuclear  quadrupole 

169 
interaction  of  Tm  in  absorbers   of  thulium  ethyl 

sulfate  (TmES).     Sources  of  ErF    (T   =  550     K;    single 

line)  were  used.     Curve  A  is  the  best  two-parameter  - 

fit  (parameters     p.     and    p    )    to  the  experimental  data. 

Curve  B  is  the  best  one-parameter-fit  (parameter     p.) 

to  the  experimental  data,    thus  disregarding  the  lattice 

contribution  to  the  electronic   shielding  (i.e.  a ,  =  y       =0), 

6  v  Z       '  CD         ' 

The  CEF    parameters  in  set    3  of  Table  II  were  used  in 

both  curves  A  and  B.     Observe  that     (AE)_     *■    0 

\       /T— *■  oo 

for  curve  B. 

The  difference  between  curves  A  and  B  shows  the 

large  contribution  of  the  lattice  part  ( 1  — v     )  q. .  to 

°  r  *co       11 

the  electric  field  gradient  at  the  nucleus.      Curve  B 
illustrates  in  particular,    that  it  is  not  possible  to  ob- 
tain a  good  fit  to  the  experimental  data  by  merely 
adjusting  the  theoretical  value  of  either    Q    or     R/~v- 
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Fig.    11:  Temperature  dependence  of  the  nuclear  quadrupole 

169 
interaction  of  Tm  in  Tm,Ov     Sources  of  ErF.. 

(T    =  550    K;     single  line)  and  absorbers  of  Tm?0, 

(5  mg/cm     of  thulium)  were  used  for  temperatures 

of  the  absorber  in  the   range    1  1°  K  <    T    <     700°  K. 

Absorbers  of  TmES  (T  =  300    K;  single  line;  5  mg/cm2 

of  thulium)  and  sources  of  Er    O     were  used  for 

temperatures  of  the  source    in  the   range    T  >  700     K. 

Curve  A  is  the  best  two-parameter -fit  (parameters 

p       and    p    )  to  the  experimental  data,    using  the  CEF 

parameters  of  Table  VI.      The  insert  shows  a  typical 

spectrum.      The  importance  of  the  lattice  contribution 

(  1  -y     )  q!Lat)  to  the  total  field  gradient  (  1  -y     )  q.(Lat)  + 
x        rco    Mn  ^  roo     Mn 

(4f) 
(  1  —  R/-0  q..         at  the  nuclear  sites  is   strikingly  demon- 
strated by  the    fact,    that  the  quadrupole  splitting    yAEL 
does  not  approach  zero  in  the  high  temperature  limit, 
but  rather  goes  through  a  minimum  and  then  increases 

again,    with       (q.)  T    — >-   ( 1  -v     )  q..  . 

0  \    11/  T — *-  00  'co       11 
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axially  symmetric  electric  field  gradient  at  the  nucleus,    i.e. 
<?3        ~  ^      ^T   =  ®'      ^n  tn*s   case  the  quadrupole  splitting      <(aE^)t  of 
the  gamma  lines,    Eq.(Z3),    reduces  to 

<AE>T     =     iC^[Pl<3^-J2)T     +4P2]  <Z4> 

Method   1:     In  order  to  obtain  the  quantities   C-,     and  <3  J  2  -  J  2  /  ^ 

Z  \  ~~z    —  /  x 

entering  in  Eq.(24)  we  use  different  sets  of  optically  determined 

3  + 
CEF  parameters  given  in  Table  II.     Set   1  was  obtained  for  Tm         in 

LaES  by  Wong  and  Richman  (14),    who  employed  observed  optical 
levels  from  a  series  of  different  optical  multiplets.     Set  Z.  was  ob- 
tained for  Tm         in  TmES  by  Gruber  and  Krupke  (15),    who  again 
used  observed  optical  levels  from  a  series  of  different  optical  multi- 
plets.    In  contrast,    set  3  was  obtained  by  a  least-squares  method 

3 
using  only  levels   observed  by  Gruber  and  Krupke  (15)    in  the      H, 

term  of  Tm         in  TmES.      The  evaluation  of  the  C       given  in  set  3 

n     6 

3 
thus  does   not  employ  optical  terms   other  than     H,    and  therefore 

should  be  the  set  most  appropriate  for  our  reduction  of  the  nuclear 
quadrupole  measurements.      To  permit  a  check  on  the  intrinsic  con- 
sistency obtainable  by  using  one  set  of  CEF  parameters  for  the  whole 
series  of  optical  levels  we  confront  in  Table  III  observed  and  cal- 
culated CEF  levels.     The  overall  agreement  is  rather  encouraging, 
the  average  deviations     between  calculated  and  observed  values  being 
only  of  the  order  of  experimental  uncertainties.     Table  IV  gives  for 
set  3  the  wave  functions  and  field  gradients  for  the  CEF  levels  which 
are  necessary  for  the  evaluation  of  the  temperature  average 
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TABLE  II 


m  -  1 

CEF  parameters  C       for  thulium  ethyl  sulfate  (units  cm      ) 

n 


c    ,  u  r°  r°  r°  <~6  d    t 

Set  number  C-  O  O,  C,  References 


1 

129.  8 

-71.  0 

-28.  6 

2 

135.2 

-71.  3 

-28.  8 

3 

130.  5 

-65.9 

-28.  6 

-•  r.    /  a  ~>  -•    o         Wont"  and 

28.6  432.8  .H  .... 

Richman  (14) 


Ciruber  and 
Krupke  (  15) 


428.  1 

427. 3         See  text 
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TABLE    III 

Observed  and  calculated  CEF  levels  for  thulium  ethyl  sulfate  in 

3  - 1 

the     H,   term  of  the  ground  multiplet  (units  cm      ).     The  calculated 

levels  of  set   1  were  taken  from  Wong  and  Richman  (14).      For  sets  2 

and  3  the  following  reduced  matrix  elements  were  used  (15): 

(j  I <y  I  j)      =   1.  0197     •    10"2  ,       (j    J 0  |  j)       =1.  5938  •    10'4      and 


(j  |y  |  j)      =   -5.  5318  '    10 


a  b 

Observed  levels  Calculated  levels 


Set   1  Set  2  Set  3 


302.  5 
274.0 


198.9 

157.  3 

110.9 

32.  1 

0 


306.  8 

304.  7 

300.  8 

281.  1 

279.  7 

274.  3 

219.  3 

221.  2 

221.  7 

212.9 

215.  1 

215.5 

204.  3 

204.  0 

198.  8 

162.  1 

161.4 

157.  8 

113.  3 

111.5 

1  10.7 

32.  1 

28.9 

51.  0 

-0.  5 

-4.4 

0.  7 

Optically  determined  levels  of  Gruber  and  Krupe  (15) 

Calculated  levels  using  the  CEF    parameters  given  in  Table  II, 
The  center  of  gravity  of  the  calculated  levels  is  adjusted  to 
give  the  best  fit. 
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TABLE  IV 

Energies,    wave  functions  and  electric  field  gradients  of  the  CEF 

3 
levels   of  the      H,   term  of  the  ground  multiplet  of  thulium  ethyl  sulfate 

(C,,    symmetry),    using  the  CEF  parameters  of  set  3  and  the   reduced 

matrix  elements  given  in  the  caption  of  Table  III. 


Energy        Degeneracy 


Wave  Function 


3J2  -  J' 
z 


137.  1 

1 

-0. 707  | 

110.6 

2 

-0.446  | 
0. 895  | 

58.0 

1 

0. 697  | 

51.  8 

1 

-0. 707  | 

35.  1 

2 

-0. 305  | 
-0. 953  | 

-5.9 

1 

0. 707  | 

-53.0 

2 

0. 895  | 
0.446  | 

131.7 

2 

0. 305  | 
0. 953  | 

163.0 

1 

0. 119  | 

-3)  +    0. 707     |+3)  -15. 0 

-2)  +    0.  895     1+4)  -1.1 

-4)  -   0.  446      I  +    2) 

-6)  -  0.  l68|o)+    0.  697     +6)         63.0 

-6)  +    0.  707   |+    o)  66.  0 

-l)  +    0.953    |+    5)  26.  3 

-5)  +    0.  305   |+    l) 

-3)  +    0.  707    |+3)  -15.0 

-2)  +    0.446    |+    4)  -ZZ.  9 

-4)  +    0.  895    |f    2) 

-5)  +    0.  953    |+    l)  -32.  3 

-l)  +0.  305    |+    5) 

-6)  +    0.986    |o)   +  0.  119I+  6)      -39.0 


a         The  general  form  of  the  wavefunction     is  given  in  Eq.    (9) 
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^3J2-J2>-,,    Eqs.    (16)  and  (24). 

A  summary  of  the   reduced  data  obtained  by  combining  our 
measurements  of  the  temperature  dependence  of  the  nuclear  quadru- 
pole  interaction  for  Tm        in  TmES  with  the  results  of  optical 
measurements  performed  on  the  same  compound  is  presented  in 
Table    V.      The  dimension^less  parameters     p       and     p       presented  in 
Table  V  are  experimentally  obtained  quantities  (compare  Fig.    10 
curve  A)  which  hold  within  the  framework  of  the  CEF  model.      The  ad- 
vantage of  introducing  these  parameters  is  that  their  deduction  does 
not  depend  on  a  knowledge  of  the  radial  distribution  of  the  4f -elec- 
trons or  the  value  of  the  nuclear  quadrupole  moment.     Such  a  know- 
ledge,   however,    enters  into  the  evaluation  of  the  shielding  factors, 
Eqs. (22),    (21)  and  (5). 

Method  2:    At  elevated  temperatures  the  temperature  average 
\3  J  2   -  J  2X,     entering  in  Eq.(24)  may  be  approximated  by  Eq.(I-17)  of 
Appendix  I,    which  yields  <^3J2- _J^      =-14.  1   C    /kT.      Expansion 
(1-17),    which  applies  to  the  case  of  an  axially  symmetric  field  gra- 
dient,   was  first  given  by  Elliott  (39).      Details  are  given  in  Appendix  I, 
which  also  includes  an  extension  to  the  case  of  non-axially  symmetric 
field  gradients.  .  From  a  plot  as  a  function  of  1 /T  of  our  measure- 
ments obtained  for  TmES  at  temperatures  T    >  200    K,    Fig.     12,    we 

obtain  from  Eqs.  (24)  and  (1-17)  the  values     P,(C^)     =  (0.  1  8  ±  0.  05)cm" 

0  -3-1 

and    P->C,   -  (2.  8±  1.  1)  10        cm      .      These  values  may  be  compared 


(38)  M.C.    OlesonandB.    Elbek,    Nuclear  Phys.   J_5,    134(1960) 

(39)  R.J.    Elliott,    Rev.    Mod.    Phys.    3_6,    385     ( 1964)  and  private 
communicat  ion 


54- 


tu 

en 

^ 

1) 

M 

<u 

M 

ti 

O 

CU 

* 

rt 

CD 

►rt 

M 

<— 1 
0 

rt  u 

CU 

•»-> 

rt 

4-> 

Q. 

rt 
•a 

INl 

■  rt 

CO 

CU 

p, 

d 

M 

X 

rt 

IM 

u 

X 

»— i 
rt 
■»-> 

Q. 

X 

X 
H 

^ 

00 

co 

in 

o 

CO 

o 

00 

o 

> 

CU 

X 

d 
rt 

d 

X 

c 

rt 

CU 

ft 

(X 

< 

1— < 

d 
cr 

^8 
i 

ro 

Tf 

m 

m 

• — < 

»—4 

Q. 

rt 

cu 

Q. 

0 

rt 

rt 

CU 

— 

00 

cu 

to 
cu 
cu 

<u 

en 

CU 

4-1 

en 

d 

0 

•  H 
4-> 

D 

4-> 

+-> 

d 

0 
X 
00 

u 
d 
d 

CU 

rt  a 

O 

CM 

O 

o 
o 

*-H 

m 

xD 

X 

■*-) 

o 

m 

rt 

X 

•M 

0 

♦J 

cu 

6 

rt 
u 

rt 

e 

v-i 

o 

l 

rt 

iTi 

mi 

in 

m 

m 
i 

•»-> 

d 

£ 

P< 

- 

rt 

o 

d 

en 

1 

i 

i 

i 

o     in 

0 

CO 

rO 

Ph 

ft 

0 

d 

en 

O 

o 

o 

o 

— ' 

i 

o 

£ 

en 

<N] 

4) 

ft 

X 

CU 

■*-> 

CU 
• — 1 

o 

<U 

r—i 

c 

cr 

X! 

4-> 

rt 

CD 

X 

£ 

CU 

d 
cu 

£ 

• 

00 

o 

. 

o 

w 

rM 

X 

QJ 

>-l 

in 

m 

in 

o 

. — . 

p. 

V) 

C 

CU 

H 

XI 

0 
43 

4-> 

00 

d 

en 

d 

CU 

d 

en 
rt 

CU 

E 

d 
o 

£ 

<u 
»-. 
d 

en 
rt 

Q. 

o 
o 

o 

o 

o 

o 

oo 
o 

o 

-H 

INI 
O 

d 

rt 

X 

rt 

£ 

^ 

+-> 

CU 

en 

i-i 

rt 

■H 

■H 

-H 

-H 

m~4 

-H 

d 

> 

3 

H 

<1 

0 

£ 

O 

d 

en 

d 

•rt 

en 

•rt 

d 

d 

cu 

£ 

'rt 
u 

00 

o 

C7* 

(M 

o 

O 

^1 

rt 
<u 
i— i 

V 

«-> 

tn 

o 

u 

•rt 

d 

w 

ro 

en 

rt 
u 

ft 

rt 

•rt 

E 

o 
u 

%-4 

P, 

o 

m 

rt 

CU 

o 

0) 

a) 

■»-> 

o 

00 

cu 

rd 

+-> 

nD 

^o 

vO 

i 

X 

(NJ 

m 

X 

rt 

4-1 

1 

1 

1 

o 

in 

■>-> 

< 

E 

H 

d 

rt 

+-> 

X 

< 

d 
•^ 

o 

o 

o 

O 

— ' 

i 

u 

H 

d 
cr 

en 
i 

c 

>. 
»— t 

u 

rt 

f— 4 
•r-l 

£ 

•iH 
CO 

r— ( 

rt 

4-J 

d 

en 

•rt 

rt 
rt 

d 

•  rt 

(0 

• 

• 

• 

1 

O 

o 

o 

0 

s-l 

d 

X 
C 

rt 

CO 

W 

6 

4-> 

en 
CU 

i—i 

rt 

X 

CU 

£ 

•rt 

u 
<L> 
ft 
X 

CU 

CU 

Ofl 
.■-1 

CU 

rt 
P. 

E 

o 
u 

en 
rt 

0 
u 
rt 

CU 

>^ 

d 
o 

CU 

•  rt 
*-• 

d 
•i-i 

rt 
•«-> 

w 

cu 
u 

d 
d 

Q. 

m 
o 

o 
-H 

in 

o 

o 
■H 

r- 

o 

in 

o 

CD 
-H 
rvj 

m 
■H 

o 

•-rt 

O 

o 
■H 

cr- 
o 

-H 

Q  =   1.5  bar 
d  Elbek  (38) 

C 

X 
O 

-1 

o 

CU 
l-l 

0) 

X 

1 

1 

«   d 

+ 

X 

4-> 

X 

c 

rt 
o 

-<-> 

rt 

d 

^    d 

£ 

H 

£ 

o 

CO 

cr 

•rt 

X 

•rt 
.—1 

rt 
> 

M 
•rt 

iTi 

X 

d 
rt 

rrt 
O 

d 

•  rt 

o 

■»-> 
<U 
CO 

- 

r\j 

m 

m 

1 

i 

>     0 
,-h    en 

rt    cu 

<->  ^ 
^  o 

O 

V-l 

rt 

a; 

(A 

rt 

en 
00 

W 

00 

d 

en 

d 

d 

o 
X 

X 

0 
X 
*-> 

»M 

(N) 

CU      r 

rt  5 
o  o 

CU     )h 

rt 
XI 

X 
CU 
CJ 
3 

CU 
•♦-> 

d 

d 

d 
cu 
> 

■r-< 
00 

•rt 

•4-> 
*J 

*rt 

>> 

X 

CU 

X 

o 

CU 

i-H 
O 
(J 

d 

•  rt 

X 

d 
rt 

d 
o 

•rt 

cu 
X 

d 

to 

CO 

CO 

co 

0 

0 

sing  the  th 
f  Tm-169  f 

X) 

X 

X 

DO 

4-> 

d 

W 

W 

W 

W 

rg 

<NJ 

d 

CO 

CU 

d 

d 
rt 

X 

(J 

rt 

o 
P. 

£ 

E 

E 

E 

E 

£ 

•c-4 

rt 

4-> 

d 

•r-l 

rt 

rt 

0) 

4-> 

cu 

£ 

H 

H 

H 

H 

H 

H 

P    0 

■*-> 

•ft 

+J 

CU 

rt 

■•-» 

o 

X 

o 

o 
ft 

Xt 

o 

■M 

4-> 

d 

■rt 

O 

rt 

-55- 

169 
Fig.    12:  Nuclear  quadrupole  interaction  of  Tm  in  absorbers 

of  thulium  ethyl  sulfate  and  thulium  oxide  plotted  as  a 

function  of  1 /T  in  the  high  temperature  ranges  where 

method  2  is  applicable  (see  text).      The  straight  lines 

are  the  best  fit  to  the  experimental  data  points. 


-56- 


absorber:  Tm2  03 
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with  the  values  obtained  by  method   1  :     p    (C    )2  =    (0.  1  57  ±  0.  00  1 )  cm  " 
and     p?C=(2.57±0.07)10        cm       .      The  agreement  between  these 
two  values  obtained  by  two  different  methods  gives  confidence  in  the 
consistency  of  our  analysis.      In  particular  we  conclude  on  this  basis, 
that  our  results  are  not  seriously  influenced  by  any  temperature  de- 
pendence of  the  CEF  parameters   C      ,    within  the  temperature   range 
studied  (9.  6    K  -  300    K).      The  agreement  obtained  for  the   results  of 
methods    1  and  2  indicates  that  the  CEF  parameter  C     is   reasonably 
independent  of  temperature.     The  higher  order  parameters  C  .,    C, 
and  C,    should  be  even  less  dependent  on  temperature,    because  of  the 
faster  convergence  of  the  associated  lattice  sums.      The  justification 
of  the  neglect  of  any  temperature  dependence  of  the  CEF    parameters 
C       in  our  analysis  is  supported  by  measurements  of  Gruber  and 
Conway  (40),    who  determined  by  optical  methods  the  energies  of  CEF 
levels  of  Tm  ions  in  TmES  at  T   =     77°  K,    194°  K  and  273°  K.      The 

changes  with  temperature  in  the  position  of  the  levels  typically  are 
less  than   10  cm      .      We  therefore  feel  justified  in  using  in  our  analysis 
one  set  of  CEF  parameters,    C      ,    determined  optically  at  a  single 


temperature. 


2.     Thulium  Oxide  (Tm.OJ 


3+  .  7X 

The  Tm  ions  in  the  Tm,0     (space  group  T,  )  occupy  two  noil' 

equivalent  lattice  sites;  sites  with  symmetry  C?  and  C,.  occur  in 


(40)         J.B.    Gruber  and  J.  G.    Conway,    J.    Chem.    Phys.    32   ,    1178 
(I960) 
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the   ratio  3:1.      Experimentally,    only  the  higher  populated  ionic   sites 

associated  with  point  group  symmetry  C?   are  observed.      By  choosing 

a  proper  coordinate  system  the  relevant  crystal  field  parameters  for 

C,   symmetry  are  limi  tedto     C,,    C,,    C  .  ,    C.    C.     ,    C  . ,    C.     ,    C ,  , 
Z     1  1  2        2        4        4        4  4        4  6 

C/  ,    C,      ,    C,  ,    C,     ,    C/ ,    C.  The  parameter  C~        has  been  elimi- 

nated by  the  proper  choice  of  the  coordinate  system.      One  may  always 

eliminate  one  C  by  such  a  choice  (22)  when  C       also  exists  . 

n  J  n 

Assume  that  in  a  given  coordinate  system  both  C       and  C  exist. 

&  7  n  n 

Then  by  rotating  the  coordinate  systems  about  the  Z-axis  by  an  angle 

6    defined  by 

tan        6  -   -  C        /C 
n  n 


one  may  reduce  the  problem  to  one  parameter  in  the  new  coordinate 
system,    C       ,    which  is  related  to  the  old  parameters  by 


c'm  =  V  (Cm)2  +  (c"m)2. 


n  n 


The  quadrupole  splitting  of  the  gamma  lines  produced  by  the  non- 
axially  symmetric  field  gradient  is   given  by  Eq.(23). 
Method  1:    Gruber  et  al.    (16)  have  studied  the  optical  absorption  and 
emission  spectra  of  Tm         ions   in   Y    O     at  the  C      symmetry  sites; 
the  same  CEF  levels  were  obtained  in  preliminary  studies  of  Tm 
in  Tm70_,    within  the  limits   of  the  experimental  accuracy.      Using 
the  energy  levels  obtained  for  diluted  Tm         by  Gruber  et  al.    (16)  we 
have  calculated  the  crystal  field  parameters  C       which  are  included 
in  Table  VI.     This  calculation  was  performed  by  minimizing      X    , 
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TABLE    VI 


3  + 


Observed  and  calculated  crystal  field  levels  for  Tm         in    thulium 


oxide  at  sites  with  C?  symmetry,    in  the     H,   term  of  the  grou 


nd 


multiplet  (units  cm      ).     The  following  set  of  CEF  parameters  was 

useda:    C°  =  -82,     C2,  =  -636,     C°A   =  -100,     C2  -  -1070,    C'Z  =  118, 
2  2  4  4  4 

C4  =  837,     C~4  =  -68,    C^  =  3,     C2  =  83,     C~2  =  2,     C4  =  227, 
C'4  =  -316,     C^  =  1,     C~6  =   154 


Calculated 
levels 

Observed 

i         i     b 
levels 

Calculated 
.Degeneracy        /3j2  _  j2\ 

Field  Gradients 

(K  +  J!> 

770 

796.9 

1                               -1.1 

-74.  2 

768 

788.5 

1                              -4.1 

-69.9 

680 

1                              16.  8 

-43.  0 

674 

1                                7.  5 

-37.9 

497 

494.  0 

1                             -    8.  7 

42.0 

429 

435.7 

1                             -21.3 

17.  5 

344 

1                             -22.2 

-    2.  3 

336 

340.  0 

1                             -17.1 

25.0 

258 

230.  3 

1                              -8.3 

11.  8 

200 

219.0 

1                                   3.  0 

22.  7 

95 

89.3 

1                             -    9.9 

42.7 

44 

30.  7 

1                               40.  3 

28.  2 

-    1 

0 

1                               25.  2 

37.5 

a        A  preliminary  set  from  Gruber  et  al.    (16) 
b        From  Gruber  et  al.   (16) 
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using  the  technique  described  by  Davidon  (4  1).     A  summary  of  the   re- 
duced data  obtained  by  combining  our  gamma   resonance  measurements 

3  + 
for  Tm         in  Tm.O^  (compare  Fig.    1  1)  with  the  results  of  the  optical 

measurements  is  included  in  Table  V. 

Method  2:     At  elevated  temperatures  the  temperature  averages 

<^3J  2  -  J^2^T     and    \J,Z,    +  ^2/>t     entering  m  Eq.  (23)  may  in  first  order 

be  approximated  by  the  expressions  given  in  Eq.(I-17)  and  (1-18)  of 

Appendix    I,    which  yield    <3J  2    -  J  2X   ='14-  1  C^/kT     and  (j  f  +   J2>~  = 

-9.4  C?/kT.      From  a  plot  as  a  function  of  1 /T  of  our  measurements 

obtained  from  Tm?0,  at  temperatures     T  =£  1270    K,    Fig.     12,    we  ob- 


tain by  using  Eqs.(23),    (1-17),    (1-18)  the  following  values: 

PlC2    [{C°i]i  +   T(C2)Z   '*    =  '°'5±0J)cm"2    and    PZ[(CZ>2+   T(C2)2 


(4. 0  ±  0. 9)  1Q        cm"  These  two  values  again  may  be  compared  with 

the  corresponding  values  obtained  by  method   1.    namely  (0.508*0.  001) 

-2  -3-1 

cm  and(3.80±0.03)10        cm       ,    respectively.      The  rather  good 

agreement  between  the  values  obtained  by  methods    1   and  2  suggests, 

as  in  the  case  of  TmES,    that  the      neglect       in  our  analysis  of  any 

temperature  dependence  of  the  CEF  parameters  C       is  a  justifiable 

approximation.     Furthermore,    the  agreement  suggests  the  absence  of 

c  r  ystallographi  c  phase  transitions   in  the  whole  temperature   range 

studied.     X-ray  diffraction  studies  of  Stecura  and  Campbell  (42)  do 

not  reveal  any  phase  transitions  within  the  temperature   range   300     K 

<T  <   1568°K. 


(41)  W.C.    Davidon,    ANL-Report  5990  rev.  ,    (Nov.    1959) 

(42)  S.    Stecura  and  W.J.    Campbell,    U.S.    Bur.    Mines,    Rept.    on 
Investigation  No.  5847  (1961) 
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VI.    ELECTRONIC   SHIELDING    FACTORS 

Our    experiments    reveal    the    presence    of    strong    charge    polari- 
zations of  closed  electron  shells.      The  shielding  (or  antishielding) 

factors  R.~,    y       and    a , ,    which  were  introduced  in  Sections  II  and  III, 
U        oo  L 

are  a  measure  of  these  charge  polarizations. 

The  antishielding  factor  y       ("lattice"  Sternheimer  factor)  may  be 
calculated  by  several  techniques  when  the  free  ion  wavefunctions  are 
known.      Wikner  and  Burns  (25),    Ray  (26),    Sternheimer  (27),    and 
Freeman  and  Watson  (28)  have  made  calculations  of  this  quantity  for 
certain  rare  earth  ions,    and  their  results  are  summarized  in  Table 
VII.      Wikner  and  Burns  used  the  (restricted)  Hartree-Fock  wave- 
functions  calculated  by  Ridley  (43)  for  Pr         and  Tm         and  calculated 
y       by  means  of  a  perturbation-variation  method.      Sternheimer  used 

the  same  wavefunctions,    but  calculated  y       by  direct  solution  of  the 

oo      ' 

inhomogeneous  Schrddinger  equation  for  the  perturbed  wavefunctions. 

Freeman  and  Watson  used  the  unrestricted  Hartree-Fock  formalism 

to  calculate  y       for  Ce       .      The  value  of  Freeman  and  Watson  for 
'oo 

Ce         is  not  very  different  from  that  obtained  by  Sternheimer  for  the 

3  + 
neighboring  ion  Pr        ,    but  differs  appreciably  from  the  value  which 

Wikner  and  Burns  obtained  for  Pr 

Theoretical  evaluations  of  the   shielding  factor  R^  are  more 
involved.      This   results  because  of  the  proximity  of  the  closed  elec- 
tron shells  to  the  distorting  source,    the  4f -electrons .      For  this 
reason,    one  may  even  expect  that  the  distorted  shells  may  produce 


(43)         E.C.    Ridley,    Proc.    Cambridge  Phil.    Soc.    56,  41  (I960) 
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repercussions  upon  the  4f-electron  shell,    as  was  pointed  out  by 

Freeman  and  Watson  (12),    (44).      Table  VII  includes  the  few  available 

theoretical  values  of  R.~  =  R       ,  +    R  for  rare  earth  ions. 

Q  rad  ang 

Theoretical  evaluations  of  the  shielding  factor    cr ?     are  physically 

similar  to  those  for  y     .      The  additional  complication  arises  from 

'co  r 

the  fact,    that  y       is  a  measure  of  the  closed  shell  distortions   ex- 
'  co 

perienced  at  the  origin,    while     cr,     is  a  measure  of  the  closed  shell 
distortions  experienced  at  the  position  of  the  4f -electrons ,    thus   re- 
quiring in  addition  a  rather  precise  knowledge  of  the  4f-electron 
density.      Theoretical  predictions  for    cr?    are  still  rather  qualitative. 
Lenander  and  Wong  (10)  came  to  the  conclusion  that  the  shielding 
factor    <r ?  was  of  the  order  of  0.  5  to  0.  75  in  the  case  of  PrCl     while 
Watson  and  Freeman  (12)  in  the  case  of  cerium  ions  likewise  con- 
cluded that  shielding  via  the    cr?     factor  is  large.      Ray  (26)  arrives  at 
the  theoretical  value  of   a ?    =  0.  52  for  the  case  of  PrCl    .      Burns  (9), 
on  the  other  hand,    using  analytic  perturbation  calculations,    con- 
cludes that  the   shielding  factor    cr      should  be  at  most  of  the  order  of 
0.  1  for  rare  earth  ions. 

The  analysis  of  our  experimental  results  yieldjthc  parameters 
p       and   p      given  in  Table  V.      Using  the  value  of  the  nuclear    quadru- 
pole  moment  Q  we  can  evaluate  the  parameters  <^  r    ^>„  and  (  1  -y     )/ 
^r2/F.     Values  of  these  parameters  are  included  in  Table  V.     It 
appeared  reasonable  to  use  a  theoretical  value  for  y       to  obtain  the 
radial  integral    ^r^^.,    since  theoretical  evaluations  of  this  quantity 


(44)         A.J.    Freeman  and  R.  E.    Watson,    Phys.    Rev.     131,    2566(1963) 
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TABLE    VII 

Theoretical  values  of  Sternheimer  shielding  factors  for  rare 
earth  ions 


R  Reference 

ang 


Ion 

y 

'  CO 

R          A 

rad 

Ce3  + 

-    73.  5 

-0.43a 

Pr3  + 

-    16.4 

Pr3  + 

-105 

Pr3  + 

-    78.  5 

T        3  + 

Tm 

-    61.  5 

1  m 

-    74 

t-,    3  + 

Eu 

(28) 
(26) 
(25) 
(27) 
(25) 
(27) 
0.29  (    5) 


a  Using  the  value  (28)     (r~     V     =4.71  a.  u. 
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appear  to  be  relatively  reliable. 

Besides  the  values  of  the   "electric"   radial  integrals    ^  r     /U  and 
^r2)^,    which  enter  in  the  analysis  of  our  measurements  of  the 
nuclear  quadrupole  interaction,    there  exists  a  "magnetic"  radial 
integral     ^r     /»#!    which  enters  in  the  analysis  of  nuclear  magnetic 
interactions.      The   effective  integral    (r      \~,  likewise  may  be  asso- 
ciated with  a  shielding  factor  (5),    (45),    which  in  analogy  with  the 
electric  case  is  defined  through  the  relation     (r      /w   -   (r      )a.(^  ~^\a^ 


(compare  Eq.(2l)).    The  difference  between  the  values  of     (r      yr    and 


(r      V     arises  because  the  contributions  from  the  closed  shells  differ 
for  the  quadrupole  and  the  magnetic  interactions.      This  difference  is 
due  to  the  different  forms  of  the  interaction  operators  for  the  nuclear 
quadrupole,    magnetic  orbital  and  magnetic  spin  interactions,    as  was 
emphasized  by  Sternheimer  (5),  (45)  and  Freeman  and  Watson  (44). 

1  1 

The  radial  integrals   (r      \~,    /r      V     and     { r2^     ,    which  enter  the 
nuclear  quadrupole,    nuclear  magnetic  and  CEF  interactions,    incor- 
porate the  contributions  to  these  interactions  from  both  the  partially 
filled  (4f)  and  the  closed  electron  shells.     These      radial  integrals  in 
principle  may  be  taken  from  experimental  observations,    a  procedure 
adopted  in  this  paper.     Table  VIII  includes  a  compilation  of  relevant 


(45)  R.    Sternehimer,    Phys.    Rev.    86_,    316(1952) 

(46)  I.    Lindgren,    Nuclear  Phys.    _3_2,    151(1962) 

(47)  B.R.    Judd  and  I.    Lindgren,    Phys.    Rev.    122,    1802(1961) 

(48)  R.  L.    Cohen,    Phys.    Rev.    134   ,    A94  (1964) 

(49)  E.    Gerdau,    W.    Krull,    L.    Mayer,    J.    Braunsfurth,    J. 
Heisenberg,    P.    Steiner,    E.    Bodenstedt,    Z.  Physik  174,  389(1963j 
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radial  integrals  for  Tm         ions  in  different  chemical  surroundings, 
obtained  from  experimental  data  by  using  the  theoretical  values   for 
Q  and  y       given  in  the  caption.      The  table  includes   radial  integrals 
evaluated  from  our  gamma  resonance  studies  as  well  as  from  other 
pertinent  experiments. 

The  interpretation  of  the   radial  integrals  in  terms  of  electronic 
shielding  factors   requires  a  knowledge  of  the  quantities      (r      ),(  and 

(  r  )     ,    as  discussed  above.      These  radial  integrals  are  not  acces- 
sible to  direct  experimental  observation  and  one  is  forced  to  use 
theoretical  values,    the  evaluation  of  which  is  presently  somewhat 
uncertain  because  of  the  lack  of  sufficiently  accurate  atomic  wave 
functions   for  rare  earth  ions.      Any  evaluation  of  electronic   shielding 
factors  is  therefore  limited  by  the  uncertainties  in  these  theoretical 
values.     Nevertheless,    by  using  a  specific   set    of  theoretical  values 
throughout  the  whole  analysis,    one  still  can  observe   the  general  trend 
in  electronic   shielding. 

Freeman  and  Watson  (50)  discuss  the  theoretical  situation  in  the 
evaluation  of        (r       )df  and     \r/df  f°r  most   rare  earth  ions.      These 
authors,    in  particular,    have   shown  that  the  values  of       (r      Vf  for 
rare  earth  ions  incorporated  in  a  solid  do  not  differ  very  much  from 
the  free  ion  values  (28).      Table  VIU  includes  a  compilation  of  elec- 
tronic  shielding  factors  for  Tm         ions  obtained  by  using  the  theore- 
tical quantities  given  in  the  caption.      The  uncertainties  of  the  theore- 
tical values  of  Q,    \T/Af  anc*   (  r      V*/-  are  presumably  less  than  30%. 


(50)         A.J.    Freeman  and  R.  E.    Watson,    Phys.    Rev.     127,    2058(1962) 
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It    appears  from  Table  VIII  that  the  shielding  of  the  nuclear  qua- 
drupole  interaction  and  the  nuclear  magnetic  interaction,    expressed 
through  the  shielding  factors  Rn  and  R*,,    is  always  small  in  the  case 


Q  -"-  "M" 


of  Tm         ions 


As  concerns  shielding  factors  other  than  R,    we  note  again  that 
our  experiments  provide  only  the  ratio  ( 1  —  y     )/   (r^      = 
(l--y     )  /  [  ( 1  -cr    )   (r2V     ],    compare  Eqs.(22),  (23).      It  appears  from 
column  7  of  Table  VIII  that  there  is  a  substantial  electronic   shielding 
associated  with  the  shielding  factor    cr?,    which  describes  the  fact  that 
the  4f  electrons  do  not  interact  with  the  direct  CEF,    but  with  a  CEF 
shielded  by  core  electrons  (primarily  5s    p     electrons).      This  ex- 
perimental observation  is  in  qualitative  agreement  with  conclusions 
we  draw  for  praseodymium  salts  from  NMR  measurements  on  lantha- 
num salts  performed  by  Edmonds  (6).      CEF  shielding  effects  of 
comparable  magnitude  were  also  obtained  by  Blok     and  Shirley  (8), 
in  the  case  of  several  rare  earth  ethyl  sulfates  and  rare  earth  double 
nitrates,    using  nuclear  alignment  techniques.     We  note  that  our    con- 
clusions concerning   <r?  are  in  agreement  with  the  theoretical 
estimates  of  Lenander  and  Wong  (10),    Ray  (11)  and  Watson  and 
Freeman  (12),    but  are  in  serious  disagreement  with  theoretical 
conclusions  of  Burns  (9). 

It  is  interesting  to  note  the  difference  in  the  a     values  presented 
in  Table  VIII  for  TmES  and  Tm    O    .      This  seems  to  indicate  that   o~  - 
depends  on  the  chemical  environment,    which  might  result  from 
different  amounts  of  overlap  of  ligand  wavefunctions  with  the  central 
rare  earth  ion.     This  seems  to  conform  with  similar  conclusions  by 
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Hutchings  and  Ray  (Zl)  for  the  case  of  PrCl,  and  PrBr    .     It  should 
be  observed  that  our  conclusions   concerning    or?  are  based  on  the 
plausible  assumption  that  y       is  much  less  dependent  on  the  chemical 
bond  than   cr ,. 

Furthermore,    we  emphasize  that  we  have  neglected  the 
non-linear  shielding  effects  (12)  in  our  analysis.     Appreciable  non- 
linear shielding  would  invalidate  the  CEF  parameterization  scheme. 
However,    due  to  the  overall  agreement  reached  in  our  analysis   --  in 
terms  of  linear  shielding  --of  the  optical  data  and  our  quadrupole 
data  we  conclude  that  non-linear  shielding  effects  play  only  a  minor 
role. 

Similar  measurements  of  the  nuclear  quadrupole  interaction  in 
TmES  as  those   reported  in  this  paper  were  performed  by  HUfner 
et  al.    (51).     HUfner  et  al.    in  the  analysis  of  their  data  did  not  take 
into  account,    that  the  optically  determined  CEF  parameter  C     does 
not  represent  an  unshielded  CEF    parameter  but  rather  represents  a 
potential  at  the  4f-electron  sites  which  undergoes  shielding  by  closed 
electron  shells  of  the  order  of  70%,    as  shown  in  this  paper.      We 
would  like  to  emphasize,    in  this   context,    that  the  lattice  contribution 
to  the  total  electric  field  gradient  at  the  nuclear  sites  is  most  easily 
observed  in  the  measurements  performed  at  higher  temperatures. 

The  nuclear  quadrupole  interaction  in  Tm,0,  has  been  investi- 
gated previously  in  a  limited  temperature   range  by  Cohen  et  al.    (23) 


(51)         S.    HUfner,    M.    Kalvius,    P.    Kienle,    W.    Wiedemann, 
H.    Eicher,    Z.    Physik  175,    416(1963) 
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and  Kalvius  et  al.    (52)    using  recoilless   resonance  absorption  of 
gamma   rays.      In  the  analyses   of  these  papers   the  direct  contribution 
from  the  lattice  to  the  electric  field  gradient,    enhanced  by  electronic 
shielding,    was  not  considered.      Preliminary  results  for  Tm.O,   in  a 
limited  temperature  range  were   reported  by  Cohen  (53). 

Although  the  importance  of  shielding  effects  expressed  by  the 
factor    cr      is     well  established,    the  absolute  values  of  the  shielding 
factor    cr      may  be  in  error  by  up  to  30%.      On  these  grounds  we  do  not 
feel  that  there  exists  any  of  the  serious  discrepancies  reported  by 
Htifner  et  al.    (51)  between  the  value  of  the  nuclear  quadrupole  mo- 
ments obtained  by  gamma  resonance  measurements  and  those  derived 
from  Coulomb  excitation  measurements. 


(52)  M.    Kalvius,    W.    Wiedemann,    R.    Koch,    P.    Kienle  and  H.  Eicher , 
Z.    Physik   170,    267(1962) 

(53)  R.  L.    Cohen,    Ph.D.    thesis,    Department  of  Physics ,    California 
Institute  of  Technology,    Pasadena,    1962,    unpublished 
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VII.    SUMMARY 
This  work  has  demonstrated  that  the  technique  of  gamma  reso- 
nance absorption  provides  a  sensitive  method  for  investigating  elec- 
tronic  shielding  by  closed  electron  shells  in  rare  earths,    via  measure 
ments  of  the  temperature  dependence  of  the  nuclear  hyperfine  inter- 
actions.     It  was   shown,    in  particular,    that  in  those  cases  where 
measurements  can  be  performed  at  elevated  temperatures  one  can 
obtain  information  on  electronic   shielding  factors  without  the  necessi- 
ty of  relying  on  CEF  parameters  determined  by  other  methods,    such 
as  optical  spectroscopy.     Our  results  lead  to  the  conclusion  that  the 
distortions  induced  in  the  closed  electron  shells  by  the  4f  shell  only 
produce     a  small  shielding  of  the  4f  electron  contribution  to  the  total 
field  gradient  at  the  nuclear  site  ("atomic"  Sternheimer  shielding 
factor        R«      <     0.  1).     On  the  other  hand  the  distortions  induced  in 
the  closed  electron  shells  by  the  CEF  lead  to  substantial  enhance- 
ment of  the  direct  electric  field  gradient  produced  by  the   surrounding 
ions  at  the  nuclear  site  ("lattice"  Sternheimer  antishielding  factor 

y     )  as  well  as  to  a  substantial  reduction  of  the  CEF  as  seen  by  the 
co 

4f-electrons  (shielding  factor  a ?) .      We  obtain  values   for  (  1  —  y      )/ 

(1—  &    )  of  250  for  Tm         ions   in  thulium  ethyl  sulfate  and  of   128  for 

3  + 
Tm         in  thulium  oxide.      The  difference  in  these  two  values   seems 

to  demonstrate  a  dependence  on  the  chemical  bond. 

It    is   interesting  to  note  that  the   ratio  of   1  -cr      for  TmES  to    1  -a 
for  Tm.O,  agrees  approximately  with  the  ratio  of  the  overall  CEF 
splittings  in  these  two  compounds. 

It  appears  that  measurements  of  the  nuclear  quadrupole  inter- 
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action  are  presently  much  better  suited  to  determine  electronic 
shielding  factors  than  to  determine  nuclear  quadrupole  moments,    due 
to  the  existing  uncertainties  in  the  different  electronic  shielding 
phenomena. 
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APPENDIX    I 

In  what  follows  we  derive  an  approximation  for  the   electric   field 
gradient  which  holds  at  elevated  temperatures. 

The  relevant  matrix  elements  entering  the  expression  for  the 
electric  field  gradient  introduced  in  the  text,    Eq.(15)  can  be  ex- 
pressed in  terms  of  spherical  harmonics  by 

(jikl|j)<3J2z-J%  =4vTW5(z    Y°(  0.,    ,.))    T  (I-la) 

i 

d-lb) 

<J  \\a\\  j)  <J*  +  J2_)  T  ^^^(sfY^.,   „.)  +  Y"Z(  ,>.,   ^.)]>T 

i 

where    the  £.   extends  over  all  4f  electrons. 
1 

Using  the  density  matrix  formalism  the  thermal  average  of 
the  spherical  harmonics  in  Eq.(I-l)  may  be  written  as 


(Ot  =  z"1  2(^MlY™exp  [-v (Ho+  v)HxM)         <j-2) 


where 


Z=     -X](xM|exp     [-)3(H     +    V)]|\M>  (1-3) 

XM 


and  /3  -   1/kT.     The  Coulomb  and  spin-orbit  interactions  are 
represented  by  H_     and    V     is  the  CEF  potential  defined  by 
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i  n  m 

and  the  electronic  wave  functions     |  X  Vy        for  pure  Russell-Saunders 
coupling  are  defined  by 

(HQ  +   V)    |  \v)       =  (Ex  +   Ey  )     |xi/>  (1-5) 

The  quantum  numbers  o-LSJ  are  represented  by  A  and  V  is  a 
quantum  number  characterizing  the  CEF  levels.  Since  the  trace  of 
a  matrix  is  invariant  to  the  choice  of  basis  functions,  we  choose 
eigenfunctions  of  J  in  Eqs.  (1-2)  and  (1-3)  rather  than  using  the 
eigenfunctions  in  Eq.(I-5)  which  are  mixed  in  M.  The  "lattice 
sums"  (2.  introduced  in  Eq.(I-4)  are  linear  functions  of  those 
used  in  the  text  (compare  Eq.(l)).      We  have  for  instance 


d  °  =  4V~^T5"  A°  (I-6a) 


a\  +CL~ZZ  =  4^2ir    I\5'a\  (I-6b) 


and  we  choose    CLc\  ~  0- 

According  to  Van  Hove  et  al.    (54)  the  exponential  factor  in 
Eqs.  (1-2)  and  (1-3)  may  be  expanded  as  follows 


(54)  L.    Van  Hove,    N.  M.    Hugenholtz,    L.  P.    Howland,    Quantum 

Theory  of  Many  Particle  Systems  (W.A.    Benjamin,    Inc., 
New  York  1961)  p.  82 


74 


CI) 

*p    [-3(H0+  V)]    =     £        pn 


d-7) 


ii  -  U 


wh 


ere    pn  =  exp  (  —  ft  H    )  and  for      n  >  0    we  have 


3  0, 

pn=(-Dny^i/ 


^ 


/  d0nexp[-(p-^l)Ho) 


0 


0 


-exp[-(i31-02)H()]  V...V  exp[-(^nl-3n_2)H0   Vexp^jy 

For  a  temperature  large  compared  with  the  CEF  interaction  energy 
(i.e.    j3E  y  <  1 )  only  the  first  few  terms  of  Eq.(I-7)  need  be  considered. 
Hence  Eq.  (1-2)  reduces  to 

d-8) 

(0TrZ''    2(xm|y^|x»m')    (x'm- |p0+Pl  +  o(/32)  |xm) 

A  A 
MM1 

where 


Z   = 


X](xm|p0+  Pl  +  o(ftz)  |am) 

AM 


and 
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(\'M»|p0|xM)       =exV(-pEx)     6XA,    6 


MM' 


-0exp(-0E     )   (\M'  |v|\M)       for     \  =  A1 


(X'M'  |pJxm)  = 


-cxp(-/3Ex,) 


exp[^(E       -EJ]-1 


E      -  E 
X'         X 


.  (x'M»  |v|xm) 

for  X    ^   A' 


Furthermore,    if  the  temperature  is  also  small  compared  with  the 

spin-orbit  splitting  (i.e.  0(  E       —  E      )  »  1,    where  X,    represents 

Xj         XQ  1 

the  first  excited  term  of  the  ground  multiplet  and     Xn     represents 
the  ground  term)  only  the  ground  term  is  appreciably  populated. 
Thus  Eq.  (1-8)  reduces  to 


(Ot-z"l2 


MM'  L 


0<xo 


M 


rm 


x0m->    <x0m-|v|x0m)    + 


d-9) 


+  2      2      <XM 

x«^x0x 


Y 


m 


X'M')     (a' 


M1 


V    X 


oM>/Ex- 


for    n  >0.     Here     Z  =  2J_  +   1     and  we  have  chosen  E      =0.      Further 

U  XQ 

more  because  of  the  Wigner -Eckart  theorem  and  the  properties  of 
the  vector  coupling  coefficients  (the  notation  of  Edmonds'  (55)  is 
used)  we  have   used  the  following  relations 


(55)  A.  R.    Edmonds,    Angular  Momentum  in  Quantum  Mechanics 

(Princeton  University  Press,    Princeton,    1957) 
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£     (\m|y™|\m)     =    (2n+   I)"2     (  X  |Yn  |  \)  •    S  (  -  1 )  ( J  M  J -M  |  J  Jnm) 

M  M 


(-1)J     M  =  (2J+  l)"2   (JMJ  -M    J  J  00) 


S      (jJ  00 1  JMJ-M)     (jMJ-m|  JJ  nm)       =6    A6      n 
.  .      \  i  f      \  I  /  nOmO 

M 


and  therefore 


l(\M  I  Ym|xM)       =0  for     n    >    0 


d-10) 


According  to  Eq.(I-l)  we  are  only  interested  in  the  cases  of 
even    n    and    m,    for  which  we  obtain  from  Eq.  (1-9) 

<Y-,»,V))T=-UJo+1)-^n-{<X0||Yn(,,     |x„>     • 


?K*D  <N,lK<"Vlk>i0+2    ,5    |>oK<"»l*'> 

i  x'  4  xQ  L 


sK^X*' IYt<->i>  I^HV1]}-^') 


-/ 


(1-11) 


In  arriving  at  Eq.(I-ll)  we  used    Eq.(I-4)  and  the  following 
relations : 
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]£      (*nM     Ym|xlM')     (\'M»     Yq|\„M) 

T"T„.     XU  nl  '        \  D   I      U       ' 


MM" 


=    [(2n+1l)(2p+l)p         <\JYnlV>       <*1Yplk> 


]T      (_l)M+]4'     (jQJ<   nm|  J0MJ'-M')     (  J'  M'  JQ  -  M  |  J  '  JQ  p  q) 
MM1 


=   [(2n+l)(2p+l)]      i(x0||Yn||v)     <X'||Y  Jx0>     (-l)m 


6  6  A(J_J'n) 

n  p       m-q  0 


where     A  (J-J'n)  -  1     if    Jn,  J'     and     n     satisfy  the  triangular  condition 
and    A     is  zero  otherwise. 

Following  Elliott  and  Stevens  (18)  we  now  make  the  following 
correspondence  between  reduced  matrix  elements 

(j|<*||j)      =  4\l  t   /5'(aLSj||2:   Y^(  0  .)    J^LSj)    Nft  y    j  (I-12a) 


(j||o||j+    l)=  -4n/t  /5'(oLSJ  ||i:    Y2(lJ.)||aLSJ  +   l)    Nfl  j    j  +   l       (I-12b) 


(j\\a  ||J+  ^=  4\/  7T    /5'(«LSJ  J  S     Y2  (^.)|oLSJ+  2)    Nfi    j    J+  2 


(I-12c) 


where 


fij    j  =  J(J+   1)  (2J  +   1)  (2J-1)  (2J  +    3) 


(I-13a) 
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"j,J+   1      3    J  (J+   0  (J+  2)  (2J+    1)  (2J+  3) 


(I-13b) 


"j,J  +  2  =  I   (J+  l)  (J+  2)  {2J+  l)  (ZJ+  3)<ZJ  +  5) 


(I-13c) 


Finally,    by  combining  Eqs.    (1-1),    (1-6),    (1-11)  and  (I-1Z)  we  obtain 
the  following  expressions: 


<J|»|J>      OlZz-£)T-A°2(r%    *<T) 


(1-14) 


<J|«|J>     <4+i2.>T  =  T   Az  <*%•<*) 


(1-15) 


where 


0(T)=  -I  (2J  +   l)'1 

b 


t<Jl«lJ>l"  nja      HCiM^OI'Oj, 


kT 


'J±  1 


J±  1 


2  |(j||q-  J  j  ±  2)  I  2  n 


J,  J±2 


'J  ±2 


(1-16) 


The  factors         12,     .     .     and        Q  in  Eq.(I-l6)  are  obtained 

from  Eqs.    (I-13b)     and  (I-I3c)  by  changing     J     to     J-l     and     J-2, 
respectively.      The  energies  E_  and    E  are  those  of  the 

center  of  gravity  of  the  terms  nearest  the  ground  term  which  have 
quantum  numbers     J±l     and    J  ±2. 
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Applying  these  results  to  the  case  of  Tm        ,    the  effect  of  the 
second  and  third  terms  of  Eq.(I-l6)  is  negligible  (     <    1%)  at  all 
temperatures  used  in  our  experiments.      Under  these  circumstances 
we  arrive  at  the  following  high  temperature  approximations  used  in 
method  2  of  the  text  (compare  Eqs.  (23)  and  (24)) 

0222-22>T  =   -i    <2J+    D'1     C°<J|H|J>     BJfJ/kT  (1-17) 

(l\  +  l2.)T  =  ~  iT  (2J+  ""'  cz  (JIfflJ>  n  j, j  AT  C"18' 


-80- 
APPENDIX    II 

Most  of  the  computations  involved  in  the  analysis  of  the  experi- 
mental data  presented  in  this  paper  were  carried  out  with  the  aid  of 
an  IBM  7094  computer.     This  section  contains  the  listings  of  the 
Fortran  IV  computer  programs  that  were  used.     The  analysis  involved 
four  major  steps. 

1)  OPTIC  ...    A  least-squares  fitting  of  the  calculated  CEF  levels, 
using  the  method  of  Section  II,    to  the  optically  observed  levels  in  terms 
of  the  CEF  parameters. 

2)  ICARME  ...    A  calculation  of  the  reduced  matrix  elements  (compare 
Eqs.(4)  and  (1-12)  in  the  intermediate  coupling  approximation  starting 
with  the  relevant  Slater  integrals  and  the  spin-orbit  coupling  para- 
meter. 

3)  QTAVE  ...    A  least-squares  fitting  of  the  calculated  nuclear  quadru- 
pole  splitting,   using  method   1  of  Section  V,    to  the  observed  splitting  in 
terms  of  the  parameters    p.     and   p    . 

4)  HITEMP  ...    A  least-squares  fitting  of  the  calculated  nuclear  quadru^ 
pole  splitting,    using  method  2  of  Section  V,    to  the  observed  splitting  at 
high  temperatures. 

Some  of  the  programs  used  are  available  in  the  SHARE  library 
and  are  not  included  here. 
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OPTIC 


The  required  card  decks  are: 


1.  MIN 

2.  READY 

3.  AIM 

4.  FIRE 

5.  DRESS 

6.  STUFF 

7.  MATMPY 

8.  FCN 

9.  THEORY    4 

10.  QSQUAR 

11.  LEIGEN 

12.  HERM 

13.  RDM 


see  SHARE  no.    980 
ZO  ANFZ013 


see  SHARE  no.    884     PK    HMEE 

see  SHARE  No.    1359    G5     XGC     0008 


Decks    1-7  were  converted  from  Fortran  II  to  Fortran  IV  and 
deck  12  was  modified  so  that  it  could  be  used  with  a  Fortran  program. 
The  listings  of  decks   8-12  follow. 
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F  C  N  ... 

"This  subroutine  is  the  link  between  theory*  and  qsquak  to  min, 
jhejvariable  metric  minimazation  routine,  share  no.  980 
input  data" 

PAR(I)=CEF  PARAMETERS  IN  THE  FOLLOWING  ORDER  C20  ,  C22  ,  C*0  ,  C<»2  , 
C4-2,C44,C4-4,C60,C62,C6-2,C64, C6-4 ,C66, C6-6  OR  C20,C40,C60, 
C66 
~NE=N0.  OF  ENERGY  LEVELS  MEASURED 
ITIME=  MAX.  XEO  TIME  ALLOWED 

ALOWER=LOWE~R~  LIMIT  OF  RELATIVE  DIFFERENCE  USED  IN  CALCULATING 
GRADIENT 

UPPER=UPPER  LIMIT  OF  RELATIVE  DIFFERENCE  USED  IN  CALCULATING 

GRADIENT 

E(  I  )=MEASURED  ENERGY  LEVELS 
R(I)=ERROR  IN  MEASURED  ENERGY  LEVELS 

V(I)=LEVEL  IDENTIFICATION       — 

DP(I)=STEP  SIZE  USED  FOR  GRADIENT 

AB=IDENTIFICATION  FOR  PARAMETERS 
SUBROUTINE  FCN(NP,G,F, PAR,MFLAG) 
COMMON/COMTHE/NFORM,P(3,40> 
COMMON /COMQSQ/NE, El  1 00 ) , S (  1 20 ) , R (  1 00 ) 
"COMMON/COMMIN/H  "  — — _____ 


DIMENSION  H(  40,40)  ,G(40)  ,  PAR  (40)  ,  V  (  1  00  )  ,  DP  (  40  )  ,  AB  {  40  ) 
DIMENSION  4H ( 50, 22  ), BH ( 50, 1 ), FP ( 40 ), I TER (40), I  TERM (40 ),OEV( 100) 
IFtMFLAG-1  II,  1,2 
'  1"READ(5,100)  NE, I T I  ME , ALOWER , UPPER 
J_00_F0RMAT(  13  /I6/2E12.5) 

CALL  ICLOCK  (INTIME)        — — : 

I  TIME= ITIMEMNTIME 


READ(5,800)(E( I ),R(T),V( I),  I=1,NE) 
800  F0RMAT(2F20.5,A6) 

RE AD ( 5 , 200 ) ( DP ( I ) ,  T=  1 , NP ) 

200  F0RMAT(F20.5) 

~READ(5,500)  AB — 

500  F0RMATU2A6) 

READ(5,900) 
900  F0RMAT(80H   HEADING 

1  ) 

WRITE(6,700)  (DP(I),  1=1, NP) 
"TOO  FORMAT(39H0STEP  SIZE  USED  IN  C ALCULAT I NG'"GRADI ENT/ ( 3H0   8E  14  .  5  )  ) 
WRITE(6,600)  AB 
600  F0RMAT(23H00RDER  OF  PARAMETERS  "  X/C20A6)  )""" 
N0DEG=NE-NP 

2  NF0RM=0 

IF(MFLAG.NE.  I)  NF0RM=1 

DO  3  I=i,NP  " 

3  P(2, I )=PAR{ I ) 
F=0.5»QSQUAR( 1,2,2) 
NF0RM=1 

DO  10  1*1 ,NP 
ITERM( I )=0 
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ITER (  I  )  =  0 
12  P(3,I)=PAR(I)»( l.O+DPf I ) ) 

FP(  I  )=0.5*QSQUAR(  1,3,2) 

RF=ABS(  (FP(I)-F)/(FP(I  )  ♦' F)T 
IFIRF.GT.  ALOWER)  GO  TO  11 
"  DP(  I  )  =  5.0»DP(  I  ) 
ITER! I  )  =  ITER( I  )  +  i 

IF(ITER(I).LT.9)  GO  TO  12 

F  P  (  I  )  =  F 

WRITE  (6,2001)~~T,DPm    "" 

2001  FORMAT (35H0GRADIENT  FOR  PARAMEJER  SET  TO  ZERO I6,E20.8) 

"GO  TO  10 
11  IF(KF.LT.   UPPER)  GO  TO  10 
" DP( I )=0.40»DP( I ) 

ITERM(  I  )=ITERM( I  )  +  l 

1  F  (I  TERM  (I  ) '.  LT79 )~  G0~T0~T2 

WRITE  (6,2002)  I,OP(I) 
"2602" FORMAT ( 35H0GRADIENT  FOR  PARAME TER'St'T  "TO  10.    16",  "E20.8)"~ 
FP( I )  =  F  +  10.0»0P( I )«PAR( I ) 

1'0"'C0NTINUE"""""' ' 

WRITE (6, 2000) (  ITERt  i),  I  TERM!  I  ),  I=1,NP) 

"^0"00~FORMAT(3OFioC5ERiVATiVE"~Movts^»PL0ST~MTNus7^rTox78rrrriH7rrai"x^^ 

WRITE(6,700) (DP( I ) ,  1=1, NP) 
DO  A  I=l,Ngr 
4  G(  I  )  =  (FP(  I  )-F)/(DP(  I  )«PAR(  I  )  ) 
"CALL  ICLOCK(LPTIME) " 

IF(LPTIME.GT. I  TIME)  GO  TO  50 

IF(MFLAG-3)  60,51,60 
50  WRITE(6, 1000) 
1000  ""FURMAT"{28H0»«n»«»«  TIME  E  XCEEDEO  $*  *$  *$  "  )" 

WRITE(6, 1001 )  (P(2,I),  1=1, NP) 
TOOl  FORMAT (32H0RESULTS  OF  FIT  UP  TO " TH I S  PO INT/3H0X= IPSE  14. 5/ 

1  (3H0   8H14.5)) 

I       WRITE (6, 1003) 

1003  FORMAT ( 13H0ERR0R  MATRIX) 
DO  6  1=1, NP 
6  WRITE(6, 1004) (H(  I  ,  J)  ,  J=1,NP) 
"100 4  "FORM  AT  (  1H01P8E1A.5/(TH08E14.5) ')"" 

_PUNC  H  1 002 , ( P ( 2,1),  1*1, NP) 

1002  FORMAT (6E 12. 5) 
DO  5  J=l,NP 


5  PUNCH  1002,  (H(J,I),  I=J,NP) 
51  PUNCH  200,(0P(I),  I=1,NP) 
"CALL  THEORY< 1,2,2,0) 

_WR  I  T E  (  6_,  900) 

ANE=NE 

DO    30    1=1, NE 

30  DEV(  I  )=SU  )-E(  If 
STAN=0.0 

DO    31     1=1, NE 

31  STAN  =  STAN+(DEV(  I  )  )»»2 
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STAN=SQRT( STAN/ANE ) 
WRITh(6,90l)(S(I),E(I),DEV(I),R(I),V(I),  1=1, NE) 

901  F0RMAT(2F20.3,F11.  1,F9.1, 15XA6)        ___      

"WRITE  (6,  700)  (OPTOi  I'=lVNP) 
ANODEG-NOOEG 
CHI2=2.0»F/ANODEG 
WR1TF(6,903)  NODEG  ,CHI2,STAN 
903  FORMATMOHONUMBER  OF  DEGREES  OF  FREEDOM=  16/ 

134H0CHI  SQUARED  /  DEGREES  OF  FREEDOM=   1PE1^.7/ 
235H0RMS  DEVI ATION"  FOR  "EACH  LEVEL=  +0R-  0PF6. 2 ) 
DO  55  1=1, NP 
DO  55  J=I,NP 
AH( I f J )=H{ I , J ) 
5  5  AH( J, I )  =  AH{ I , J ) 

CALL  MATINVI AH, NP,BH,0, DELTA) 
WRIT£(6,  1005)"  DELTA 
1005  FORMAT  WHODELTA  =  E20.8) 

~  IF  (LPT  I  ME.  GT.  IT  I  MET"  CALL  EXIT 
60  RETURN 
END 
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THEORY  4  ... 

THIS  SUBROUTINE  CALCULATES  THE  ELECTRONIC  ENERGY  LEVELS  OF  A  RARE 
EARTH  ION  IN  NTERM  TERMS  OF  THE  OPTICAL  SPECTRA.   THE  MAX  J=l5/2. 
"HIGHER  J  MAY  BE  HANDLED  IF  THE  APPROPRIATE  DIMENSION!  STATEMENTS 
ARE  CHANGED.   THE  LEVELS  ARE  CALCULATED  FROM  GIVEN  CEF  PARAMETERS. 
NTERM  MAX  =  20,  MAX  NO.  OF  LEVELS  =  100,  ANY  POINT  SYMMETRY  THAT 
USES  ONLY  CEF  PARAMETERS  WITH  EVEN  N  AND  M. 
INPUT  DATA 

NTERM=NO.  OF  OPTICAL  LSJ  TERMS 

NJ=NO.  OF  DIFFERENT"  J  VALUES  IN^TERMS" 
XJ=MIN.  J  VALUE  -  1.0 
"IW=DUMMY  VARIABLE 

IDEG=1  FOR  INTEGRAL  J,  2  FOR  HALF  INTEGRAL  J 
"DELTA=LIMIT  FOR  OFF-DIAGONAL  ELEMENTS  IN" E IGENVALUE  SUB-" 

ROUTINE  EIGENH 

NRDM=STARTING  POINT  FOR" RANDOM'  NUMBER  "TUNCTI ON  RDM ~ 
SYMTRY=C2  OR  C3H  AS  THE  CASE  MAY  BE 
"  C3HSYM=C3H  IF  THIS  IS  THE  CASE 

THETA(N,M)=REDUCED  MATRIX  ELEMENTS,  N  IS  ORDER, M  IS  TERM 
"IDENTIFICATION 

AJ(I)=J  VALUE  FOR  THE  I TH  TERM  

NOBSVR(M,Nr=l  IF  THE'NTH  LEVEL"  IS  NOT"  OBSERVED  AND"  ZERO"" 
OTHERWISE 
"NORMAL(M) =N0.  OF  THE  LEVEL  IN  THEMTH  TERM  WHICH  IS  TO  BE" 
NORMALIZED  TO  ZERO 

NOP=NO.  OF  CARDS  WITH  OPERATOR  EQU I VALENT  MATR I  X  E L EMENTS  FOR 
EACH  J  VALUE  DIVIDED  BY  3 


HEADIK, I )=LABELING  FOR  E I GENFUNC T I QNS 
FMT(K,J)=FORMAT  FOR  E I GENFUNCT I ONS 
"OA,  OB,  OC(K,I,J)=IJ  OPERATOR  EQUIVALENT  MATRIX  ELEMENT  FOR 
KTH  J  VALUE  FOR  ALPHA,  BETA,  GAMMA  RESPECTIVELY,  I  AND  J 
GO  FROM  1  TO  2J+1,  I  IS  INDEX  FOR  MJ=-J  AND  2J+1  IS  INDEX 

FOR  MJ=J  

SUBROUTINE  THEORY  { I Q , JQ , KQ, IOP )" 

DIMENSION  0A(8,16,16),0B(8,16,16),0C(8,16,16),THETA(3,20),W(16) 
DIMENSION  CEF(16,16,2),GIVE(300),ANM(3,4,2),AJ(20),SIGN(2) 
DIMENSION  E(100),R(100),V(100),P(3,<tO),X(20),QSQ(3),DP(20),S(123) 
"DIMENSION  DEF( 16, 16,2) ,W1( 16) 

_D I  MENS  ION  HE  AD (8, 33) ,FMT(8, 12) ,FMTA(  12),EVEC1 ( 16) ,EVECR(  16) 

DIMENSION  C( 16, 16) 

DIMENSION  NOBS VR( 20, 20), NORMAL (20) , Sit  120) 
COMMON /COMTHE/NFORM,P 
COMMON/COMQSQ/NE,E,S,R 
"C0MM0N/C0MWAV/C1 ( 256) 

^COMMON  DUMMY(281),  UP 

COMPLEX  C,CTEMP,C1 
COMPLEX  CS0RT,CONJG 
CALL  0VERFL(J123) 
10  IF(NFORM)  1,  1,2 
1  REAO   (5, 100)NTERM,NJ,XJ, IW, I DEG, DELTA , NRDM 
100  FORM AT ( 21 2, F 10. 5,2  I 2/E 20. 5, 13) 
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NRDM=NRDM* I 

DO  91  I=1,NRDM 
VI  IRDM=RDM(XRDM)  

READ(5,73)  S YMTRY , C3HS YM 
73  FORMAT ( 12A6) 

DO  999  M=1,NTERM 

READ   (5,101)(THETA(N,M),N=1,3),AJ(M) 
101  F0RMAT(3E20.5,F10.5) 

RE AD (5, 99)  (NOBSVR(M,N),  N=l, 16) , NORMAL (M) 
~99" FORMAT { 161 1, 12) 
999  CONTINUE 

DO  1000  1=1,8 

DC  1000  J=l,16 

DO  1000  K=l  , 16 

OA(I,J,K)=0.  

08 (  I  ,J,K)=0.' 

0C(  I  ,  J,K)=0. 
.000  CONTINUE 


DO  998  K=1,NJ 

RE  AD  (5, 103)    NOP,  (HEAD(K,  I  )  ,1  =  1,33)  ,  (  FMT  CK,~J  )  ,  J  =  l ,  12  ) 
J03_F0RMAT(  I2/20A4/13A4/12A6) 

DO  997  L=l,  NOP '"" 

READ(5, 102)  I, J, (OA(K, I , J) ) , I , J , ( OB ( K , I , J ) ) , I , J , { UC ( K , I , J ) ) 
102  FORMAT  (212,  E20.8) 
997  CONTINUE 

AK  =  K 
_LMAX=AK+XJ+0.52 

DO  996 L=1,LMAX 

MIN=L+1 


JMAX=2.0«(AK+XJ)+1.01 
MAX=JMAX+1-L 
DO  996  N=MIN,MAX 
_NA=JMAX+1-N 
— OA  (  K  ,MAX  ,  N  )  =  0A  (  K",  NA  ,  L  ) 
OB(K,MAX,N)=OB(K,NA,L) 
OC(K,MAX,N)=OC(K,NA,L) 
996  CONTINUE 

DO  995  L=1,JMAX 
DO  850  N=1,L 
OA(K,N,L)=OA(K,L,N) 
OB(K,N,L)=OB(K,L,N) 
OC  (K,N,L )  =  OC(K,L,N) 


850 

CONT INUE 

995 

CONTINUE 

998 

CONTINUE 

c 

SETTING 

DO  1001  1=1,3 

DO  1001  J=l,4 

DO  1001  K=l,2 

ANM( I , J,K)=0. 

1001 

CONTINUE 

"UP  OF'ANM 
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2  IA  =  0 
I  S=  I 

_1SA=1 _____ 

ICTR=0 

DO  994  N=l,3 
MAX=N+1 

DO  994  M=1,MAX 
"DO  994  L=l,2 

IF(L-N)  111,111,994 

"lTl~CONTINUf  ' 

IF(M-L )994,993,993 
"993  ICTR=ICTR+l 

IFISYMTRY.EQ.C3HSYM) 
1  GO  TO  (5,994,5,994,994,994,994,5,994,994,994,994,5,994) ,  TCTR 

5  IA=IA+l  

IF( I  A- IQ) 992,99  T,  992 
992  ANM{N,M,L )=P(KQ, IA) 

GO  TO  994 
991  ANM(N,M,L)=P( JQ, IA) 
"994~*C0NTINUt 

IF(NFORM)  112,113, 112 


l'lV~WR  I  T  E  ("6",  3000  )  (  (  (  ANM  ( N",M,  L")  rT^TT^Jf  M=1,4J,  N*lf3) 
3000  FORMAT  (5H0ANML/(  6E20.8)  )_ 
~~112"CONTINUE 
C  CALCULATE  CEF  MATRIX  ELEMENTS 

SIGNt 1 ) =1.0 

SIGN(2)=-i.O     

DO  900  T=i,NTERM 
IJ=AJ{  I  )-XJ+.01 


"JMAX  =  2.0«AJ(  I  )  ♦  l'.OV" 
IFISYMTRY.NE.C3HSYM)     IA=NP-NTERM 

"  I  F  ( I  0.  GT  .  I A+  1 )   ooro"52~~ 

DO    903    J=1,JMAX 

D0"903"  K=~V,~JMAX 

OG  903  L=l,2 
' C£F( J,K,L)=0.0 
903  CONTINUE 

"DO  901  J=i,JMAX 

DO  901  K=1,J       

IF (J-K) 8 00,801,800 

801  M=l 

GO  TO  810 
800  IF( J-K-2)901,802,803 

802  M=2 

GO  TO  810 
"80  3~T  F  (  J-K-4  1901,804,805 

804  M=3 

GO  TO  810 

805  1F( J-K-6)901,806,901 

806  M=4 

810  CONTINUE 
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DO  899  L=l,2 

CEF(J,K, L)=SIGN  (DM  THETA(l,I)»OA(IJ, J, K)«ANM(1,M,L) ♦THETA(2,I )• 
1     OB(IJ,J,K)»ANM(2,M,L)+THETA(3,I)«OC( I  J, J,K)»ANM( 3,M,L) ) 

CEF(K,J,L)=SIGN(L)»CEF(J,K,L) 
899  CONTINUE 
90T  CONTINUE 

DIAGONALIZE  THE  HERMITIAN  MATRIX  CEF 

IF( IOEG-2)500,501,500 
501  JMAX=JMAX/IDEG 
~~JAB=1 

00  503  J=l,JMAX 

KAB=1 

DU  504  K=1,J 

00  505  L=l,2 

DtF( J,K,L)=CEF{ JAB,KAB,L  ) 

def  {  k  ,j  ,  lt=  siGNrrn  •DTfFrjyir,  n 

505  CONTINUE 

"KAB=KAB+2 
504  CONTINUE 

JAB=JAB+2 
503  CONTINUE 

T  F  (  J  MAX .  EQ.T)~W  (  1")  =OEF  (  I  ,T,T) 
IF( JMAX.EO.l )  GO  TO  506 
'"  IF(  JMAX.NE.2)  GO  TO  499 

TEMP  =  SQRT((DEF(1,1,1)-DEF(2,2,1))*»2  4-4.0»(DEF(2,1,1)*»2  +  DEF{2,1,2) 
1**2)  ) 

_W( 1)  =  ( DEF (1,1,1 )+DEF( 2,2,1 )+TE MP)/ 2.0     

"W  (  2  )  =  (  DE F  (  1 ,  lVl )  ♦DEF  (  2  ,  2  ,  1  )  -TE  MP  )  /  2 .  0 
GO  TO  506 

499  CALL  EIGENH(D£F,W,JMAX, DELTA) 
GO  TO  506 

500  MZ  =  0 

DO  515  M=l,2 

JMAX1=(JMAX-1 )/2 

IFIM.F0.2)     JMAX1={ JMAX*1 )/2 

JZ=2 

IFCM.EQ.2)     JZ=1 

DO    507    J=1,JMAX1 

KZ  =  2 

"   IFIM.E0.2)    KZ=1~ 

DO    508    K=l,J 

DO    509    L=l,2 

DEF( J,K,L)=CEF( JZ,KZ,L) 
509    DEF(K,J,L)=SIGN(L)*DEF(J,K,L) 

508    KZ=KZ  +  2 

507    JZ=JZ+2 

IF( JMAX1.EQ.  1)     Wl(  l)=DEF(  1,  1,  1  ) 

IF( JMAX1.EQ. 1)     GO    TO    511 

IF ( JMAX1.NE.2)     GO    TO    512 

TEMP=SQRT( (DEF ( 1,1,1 ) -DEF (2,2, 1 ) ) • *2+4. 0« ( DEF ( 2 , 1 , 1)*»2*DEF(2, 1,2) 
1**2)  ) 
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512 
51  1 

510 
515 
506 


18 


17 


69 


50 


52 
53 


51 


6005 
6006 


Wl  (1  )  = 
Wl  (2)  = 
GO  TO 
CALL  E 
DO  5  10 
MZ=MZ* 
W(MZ)= 
CONTIN 
KMAX=J 
DO  17" 
KP1=K+ 

00  17 
IF (W(K 
T£MP=W 
W( K) =W 
W(  J)  = 
CONTIN 
INORM= 
TEMP=W 
DO  69 
W(K ) =W 
DO  50 
IF (NOD 

"S  (  I  S  )  = 
SKIS) 

1  S  =  I  s  ♦ 
CONTIN 
SO  TO 
IF { IDE 
IZS=KQ 
I2T=IA 

~IF(  10. 
DO  5  1 

"IF  (NOB 
S( ISA) 
I  SA= I S 
CONTIN 
IF{ IOP 
IF( IDE 

"CALL  E 
DO  600 
DO  600 
INDEX= 
C(K, J) 
KMAX=J 
MULTIP 
DO  400 
TtMP=R 
IFITEM 
CTEMP= 
GO  TO 


(0 
(D 
51 
IG 
J 
1 

Wl 
OE 
MA 
K  = 
1 

J  = 
)- 
(K 
(J 
TE 
OE 
NO 
(  I 
K  = 
(K 
K  = 

sv 
w 
=  s 

1 

OE 
53 

G. 


EF( 1,1,1 )+DEF(2,2, l)*TEMP)/2.0 
EF(1,1, 1)+DEF(2,2, l)-TEMP)/2.0 

I  

ENH ( DEF,W IV JMAXl, DELTA) 
■If JMAX1 


(J) 


x-i 


iiKMAX" 


KP1, JMAX 

W(J)  )17, 17,  18 

) 

)     

MP 


KMAL( I ) 
NORM) 
1, JMAX 

)-t£mp 

lTiMAX  ' 

R(  I,K) .EQ.l) 

K) 

(  IS) 


GO  TO  50 


EQ.2)  JMAX=JMAX/2 


+  1 

EQ. 

K=l 

SVR 

=  S1 

A+l 

Ot 

.NE 

G.E 

IGE 

5  J 

5  K 

J  +  K 

=  C1 

MAX 

LY 

0  J 

EAL 

P.E 

CSQ 

400 


IZT)  IZS=JQ 

,  JMAX 

(  I,K).EQ.T)"GO  TO 

(  ISA)*P{  IZS, IZT) 


'51" 


.0) 
Q.2 
NH( 
=  1, 
=  1, 
♦  JM 
(  IN 
-I 
C  B 
=  1, 
(C( 
Q.O 
RT( 
2 


GO  TO  900 
)  JMAX=2*JMAX 
CEF  ,  W,  JMAX,  DELTA)" 
JMAX 
JMAX 
AX-JMAX 
DEX) 


Y  PROPER~PHASE"FACT0R~ 

JMAX 

1,  J)  ) 

.0)  GO  TO  4001 

CONJG(C( JMAX, J) )/C( 1, J) ) 
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400  1  CTEMP=CSQRT(-CONJG(C(KMAXfJ))/C(2tJ)) 

4002  DO  4003  K=1,JMAX 

4003  C(K, J)=CTFMP»C IK, J)  _ 

"4000  CONTINUE 

DO  6017  K=l,KMAX 
~KPi=K+i 
DO  6017  J=KP1,JMAX 
IF<W<K).LE.W( J) )  GO  TO  6017 

__TEMP=W(K)       

_W(K)=W( J) 
W( J)=TFMP 
DO    6018    L=1,JMAX 
CTEMP=C(L,K) 
~C(L,K)=C(L,J) 

6018    C(L?J)=CTEMP 

E01T~C0NTINUF~ 


"6000 


'6001 


6003 

6002 

900 


WRITE (6, 6000) (HEAD(IJ,K),K=1,33) 
FORMAT (33A4) 
DU  6001  K=i,  12 
FMTA(K)  =  FMT(  IJ,K)~ 

DO  6002  J=1,JMAX 

DO  6003  K"=TtJMAX 

EVECR(K)=REAL(C(K,J) ) 

EVECI (K)=AIMAG(C(K, J) ) 

WRITE(6,FMTA)  W(J)  ,  (EVECR(K) ,K= I » JMAX ) , ( EVEC I ( K ) ,K=1, JMAX) 

CONTINUE 

RETURN 

END T 
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C      QSQUAR... 

"C     'FUNCTION  TO  CALCULATE  VALUE  OF  Q»»2,  WHERE  CHI»*2  IS  MINIMUM  VALUE 
C UF  THIS  FUNCTION  

FUNCTION  QSQUAR( IQ,JQ,KQ) 

DIMENSION  E(  100)  ,R(  100),  V{  100)  ,P(3,20)  ,X(20)  ,  QSQ  (  3  )  ,  DP  (  20  )  , S( 1500) 

DIMENSION  GI VE(300) 

COMMON  GIVE,S,E,R,V,P,X,QSQ,NE,NP, IQ,JQ,KQ,IOP,XTEST 

OSOUAR  =  0.0 

CALL  THEORYt IQ,JQ,KQ, 1 ) 


DO  10  L=T,NE     ; 
10  QSQUAR=QSQUAR+t ( E ( L )-S ( L ) ) /R( L ) ) »»2 

RETURN "" 
END 
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LEIGEN... 

THIS  IS  A  LINK  BETWEEN  A  PROGRAM  CALLING  EIGENH  AND  HERM 

SUBROUTINE  E I GENH ( CEF , W , JMAX , DEL T A )   

COMMON /COMHER/b( 16) , H{200) 

COMMON/COMWAV/C 

COMPLEX  H,C,CMPLX 

DIMENSION  C(16,l6),CEF(16,i6t2),W(l6J 

IREGIN=16-JMAX 

DO    1     I=1,JMAX  

~ISM*I  BEGIN 
B(  IS)=CEF(  1,1,1) 
1  =  1 

IMAX=JMAX-1 
DO  2  J=i , 1MAX 

JP1=J+1  

"DO    2    K=JP1,JMAX~ 
H( I )=CMPLX(CEF (J,K,1),CEF(J,K,2)) 

KZE=I+I/2-l 

H(KZE)=(0. 0,0.0) 

H(KZE+1)=(0. 0,0.0)      •  

"C  A  LI  HERM(H,JMA  XTCTO  YD  E  L  T  A  ,  IT) 
DO  3  1=1, JMAX 
IS=I+IBEGIN 
W(  I  )=B( IS) 
RETURN 
END 
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•  HERM... 

•  THIS  SUBROUTINE  DIAGONALIZES  A  HERMITIAN  MATRIX  AND  OBTAINS  ALL 

•  EIGENVALUES  AND  EIGENVECTORS.   IT  IS  SHAKE  NO.  884  PK  HMtE  WHICH 
V~  ""HAS  BEEN  MODIFIED  SO  THAT  IT  MAY  BE  CALLED  FROM  A  FORTRAN  PROGRAM" 

•  ONLY  THE  MODIFICATIONS  ARE  SHOWN  HERE. 

»     '  THE  CALL  STATEMENT  TO  tiE    USED  IS  CALL  HERM { H , N , U , PR , DEL T A ,  I  T  ) 

•  ALL  ARGUMENTS  HAVE  SAME  MEANING  AS  DESCRIBED  IN  SHARE  NO.  884 

•  "WRITE  UP,   PR  SHOULD  ALWAYS  BE  ZERO  AND  IT  IS  THE  NUMBER  OF 

•  ITERATIONS 

HERM    SAVE*  1,2 

SXA  XR4.4 
"CAL  3,4 

STA  HI 

CAL»  4,4~ 

LGL  18  

ST0'~H  1 

CAL  5,4 
"STA  H2~ 
_CAL»  6,4 

LGL"18 

STD  H2 


CAL»  7,4 

SLW  H3 

TSX  HMEE,4 
HI  PZE  *♦,,»* 
H2  PZE  **,  ,"••" 
H3  BSS  1 


XR4  AXT  »»,4 

CAL  COMMON+20 
SLW«  8,4 
RETURN  HERM 
FIRST  CARD  OF  SHARE" NO.  884  FOLLOWS  THIS  CARD 


•      LAST  CARD  OF  SHARE  NO.  884  PRECEDES  THIS  CARD 
COMMON  BSS  23" 
END 
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ICARME 
The  required  card  decks  are: 

1.  ICA 

2.  RME 

3.  SlXj 

4.  DELTA 

5.  FACT 

The  programs  in  decks   3  and  4  were  written  by  B.  A.   Zimmerman, 
The  listings  of  decks    1-5  follow. 
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— row 


i 


I  C  A  ..  . 
THIS  IS  A 
IN    THE     IN 

"FOR  TM  IV 
COUPLING 
REF.  F.H. 
J.  CHEM. 
INPU 
El, 
F2, 
ZETA 

"DIMENSION 
DIMENSION 
DIMENSION 
READ(5,  10 
FORMAT  ( 
F2=( El  +  14 
F4=(6.0»F 
F6=( 7.0*F 

7ETA=-ZET 
WRITEJ  6,  1 

"FORMAT  (AH 
4H0F2=E14 
PURE^ELEC 
E3H  =  0.0 
E3F""~=~15Y 
E1G_= 
E1D  = 
EH  = 
E3P  = 
EIS  = 
COULOMBIC 
S  (  1  ,  1  )  = 
S(2,2)  = 
S(3,3)  = 
S(2,  1)  = 
S(3,  1)  = 
S  (  3 ,  2  )  = 
D0_2  I  = 
DO  2  J  = 
S(  J,  I  )  = 
A  S  (  1  )  =  1  .  0 
AS(2)=0.0 
"AS  (3)*  1.0' 
AL( 1 )  =  3.0 
AL (2 )=2.0 
AL(3)=1.0 
AJ=2.0 
NORDER  = 
I  FLAG  =1 
WRITE  (6, 


-  5 
~44. 

50. 
70. 
85. 


PROG 
TERME 

OR  P 
PARAM 

SPED 
PHY. 
T  DAT 
E2,  E 
F4T  F 
1  =  SPI 
"AS  (3 

S(3, 

ERAS 
0)  £1 
F20.5 
3.0»E 
2-39. 
2-18. 
Al/2. 
000) 
1  E 1  =  1 
.7,5X 
TROST 

0  •  F 
.0  • 
0  «  F 
0  *  F 
0  ♦  F 
0  •  F 

AND 
E3F  ♦ 
E1D~ 
E3P  - 
2.0  • 
0.0 
-3.0 
2 ,  3_ 
If  2 
S(It  J 


RA 

DI 

R 

ET 

DI 

32 

A 

3  = 

6 ~ 

N- 

r, 

3) 

S 

,E 

T~ 

2  + 

0« 

0* 

0 

El 

PE 

3H 

AT 


M  TO  CALCULATE  THE 
ATE  COUPLING  APPRO 
IV  GIVEN  THE 
ER,  AND  PURE 
NG,  PHY. REV. 
,  1178(1960) 


OPERATOR  EQUIVALENT  FACTORS 
XIMATION 
SLATER  INTEGALS,'  SPIN-ORBIT"- 
L-S  OE  FACTORS. 
58,255(1940)  AND  GRUBER  AND  CONWAY 


RACAH'S  LINEAR  COMBINATION  OF  SLATER  INTEGRALS 


ORBIT  COUPLING  PARAMETER 

AL(3) 

,  EV(3,3),  E ( 3  ) 

(20) ,A(3),B(3),G(3),AE(3),BE(3),GE(3) 

2,E3,ZETA1 


li.0»E3)/42.0 

E2-E3)/ll.O 

E2-3.0«E3)/77.O 


,E2,E3,ZETA1,F2,F4 
14.7,5X3HE2=E14.7, 
F4=E14.7,5X3HF6=E1 
I C~  ENERGY  LEVELS ~ 


2  V "18.0  •  F4  -  273. 
F2__+  14  8.0  »  f_4  ♦  91 
2  -  48.0  *  F4  ♦  72  8. 
2  +  60.0  •  F4  ♦  14.'0 
2  +  84.0  '•  F4  -  12  74 
2  «■  249.0  «  F4  +  172 
SPIN-ORBIT  MATRICES, 
4.0  »  ZETA 


,F6,ZETA 

5X3HE3  =  E  1 47775  X  6HZTTA  »^E"  FSTTT" 

4.7, 5X5HZETA=E14.7) 


0  •  F6 

.0  •  F6 

0  *  F6 

»  F6 
.0  •  F6 
9.0  »  F6 

Characterized  by" j 


ZETA 

2.449490  •  ZETA 


•  1.414214  •  ZETA 


3 

400) 
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400  FORMAT  (  1H1  12X  6HENfcKGY  15X  3H3F2  15X  3H1D2  15X  3H3P2  13X 
15HALPHA  14X  4HBETA  13X  5HGAMMA   ) 

GO  TO  30 

1 0  ' S (  1 •  1 )  =  E3H  ♦  6.0  »  ZETA 

S(2,2)  =  EIG 

~S(3,3)  =  E3F  -  3.0  •  ZETA 

S(2,l)  =  2.0  »  3.162278  /  1.732051  •    ZETA 

S( 3,  1)  =  0.0 

S(3,2)  =  -  2.0  »  3.316625  /  1.732051  «  ZETA 

S(l,2)  '"'=  S(27H 

S( 1, 3)  =  0.0 


S(2,3)  =  S(3,2) 
AL( 1 )=5.0 
AL(2)=4.0 

AL  (  3  )  =  3.0        

AJ=4.0 
NORDER  =  3 
IFLAG  =  2 
WRITE  (6,500) 
5~00  FORMAT  (  1H0  12X6HENERGY  15X  3H3H4  15X3H1GA  15X  3H3FV 13X 

15HALPHA  14X  AHBETA  13X  5HGAMMA   ) 

GO  TO  30~ 
20  S( 1, 1)  =  Ell 


S(2,2)  =  E3H  -  5.0  *  ZETA 
S(2, 1)  =  -  2.449490  •  ZETA 
S  (  1 ,  2  )  =  S  (  2  , 1  ) 

_AS( 1 )=0.0 

AS(2)=f.~0 
AL ( I )=6.0 
~AL(2 )=5.0 
AJ=6.0 
NORDER  =  2 

IFLAG  =  3      , 

WRITE  (6,600) 
600  FORMAT  ( 1H0  12X  6HENERGY  15X  3H1I6  15X  3H3H6  13X 
"  15HALPHA  14X  4HBETA  13X  5HGAMMA    ) 

TEMP  =  SQRT((S(1,1)+S(2,2))*»2*4.0MS(1,2) )«»2-4.0»S( l»l)»S(2,2)) 
E(l)=(S(l,l)+S(2,2)+TEMP)/2.0 

_E12)=(S( 1, iH-S(2,2)-TEMP)/2.0  

DO  12  1=1,2 

TEMP=(E(  I )-S(  1,  1  )  )/S(  1  ,2  ) 
EV(  1,1  )  =  i.O/SORT(  1.0*TEMP»*2) 
12  EV(2, I )=TEMP/SQRT( 1.0+TEMP»»2) 
GO  TO  31 

30  CALL  E IGVV(S,EV,E, NORDER, ERASS)      

3100  25  1=1, NORDER 

A( I )=RME( AS( I ) ,AL( I ) ,AL( I ) , A J, A J, 2.0) 
D(  1  )=RME(AS(  I )  ,AL( I  ) ,AL(  I  )  ,AJ, A J, 4.0) 
25  G(I)=RME(AS(I),AL(I),AL(I),AJ,AJ,6.0) 
1F(N0RDER.EQ.2)  GO  TO  5 
A13=RME(AS(1) ,AL(1), AL ( 3 ) , AJ , A J , 2 . 0 ) 
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AS( 1 ) ,AL( L ),AL (3) , A J, A  J, 4.0) 
ASH)  ,AL(1)  ,AL(3)  , A J, A J, 6.0) 


)  (A(  I  )  ,B(  I  )  ,G(  I  ),     I=1,N0RDER),A13,613,G13 
8) 

NORDER "' 


'NOR0£R 
I )     ♦    A( J) 


(EV( J, I ) )»»2 
( EV( J, I  )  )»»2 


I  )     ♦    B( J) 

I )     ♦    G( J)     •     (EVU,  I  )  )  *«2 

Q.2)  "GO    TO    3 

I  )     ♦    A 1 3 

rrvBiT 


I  )    ♦    Gi3 


EV(  I,  I  ) 

EV(  LVIT 

EV(1,I  ) 


EV(3,  I )»2.0 

El/T3,n»2.0' 

EV(3,I )»2.0 


R.EQ.2)     GO    TO    63 

,200)  (E(  I  )','  (EV(  J  ,"I  )7J    =""1,N0R0ERT, 

ORDER)  

LH0i>7eT8TTl 

0,20),     IFLAG 

,201  )  (E(  I  ),     (EV{  J, "I  )Vj"=  "lVNORDER)  , 

ORDER) 

1H01P6E18.7T 


AEU)V  BETC  I),  ~0E[l  ), 


AEU  )V  bE(l) ,    GL(I ), 
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R  M  E  . 
THIS  IS 
OPERATO 
"THE  ELE 
REE.  B. 
THIS  VE 
FUNCTI 
"  TEMPJ= 
TEMPL= 
" L1=AL~ 
L2  =  BL 
J1  =  AJ 
J2  =  BJ 
IS  =  AS 
WI=SIX 
"" W2=SIX 
N  =  AN 
'"  J  =  AJ 
SFACT= 
"~RME1=( 

GO  TO 
2~RME=16 
RETURN 
4  ~RME=-3 
RETURN 
6~RME=12 
RETURN 
i~  WRITE"  (' 
0  FORMAT 
STOP 
END 


A  FUNCTION  TO  CALCULATE  THE  REDUCED  MATRIX  ELEMENTS  OR 
R  EQUIVALENT  FACTORS  FOR  THE  4F12  ELECTRON  CONFIGURATION. 
MENTS  ARE  OF  THE  FORM    (S,L,J  11N11  S , L  '  V  J  '  )  . 
R.JUDD,  PROC.  ROY.  SOC.  A241 , 4 14 { 1957 ) . 
RSION  IS  FOR  J=J*. 
ON  RMEI AS,AL,BL,AJ,BJ, AN) 
SQRTI  (2.0»AJ-H.0)*(2.0«BJ+1.0)) 
SQRTI (2.0«AL*1.0)»(2.0*BL*1.0)) 


J(AL,AJ,BL,BJ,AS,AN)»(-1.0)»»(L1*-L2«-JI*J2) 

J  I  3 . 0 , AL , 3 . 0 , BL , 3. 0 ,  AN  )*(-!. 0)**<L  1>1 2VJ  1 «TJ7T 


S0RT(FACT(2*J-N)/FACT(2»J  +  N+1)  ) 

-1.0)»* ( IS-  J +2   )*TEMPJ»TEMPL»SFACT»W1»W2 

(  1,2,1,4, 1,6) ,N        

.~0»SQRT(7.0/15.G)*RME1 


2.0*SORT( 14.0/11.0)»RME1 
80.0»S0RT(7.0/429.0)»RME1 

6, 106) 

(32HORME  DOES  NOT  EXIST  FOR  THIS  N 
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C      SIXJ.. . 

C  A  FUNCTION  FOR  SIX  J  SYMBOLS        2/25/64 

C ROTENBERG  ET.  AL .  PAGE  13  EOUATION  (2.3) 

~C_        'INPUT  SIXJ( Jl, J2,L2,L1,J3,L3)  IN  FLOATINGPOINT 
C       REQUIRES  DELTA  AND  FACTORIAL  ROUTINES 
"FUNCTION  SIXJ(AfB,C,D, E.F) 
TRI1  =  A+B-E 
IF(TRH)2,l,r 

2   SIXJ  =  0.0    

RETURN 
1   TRI  I  =  A-B+E 

IF (TRI  1)2, 3, 3 
3   TRI  1  =  -A-H3  +  E 

IF (TRI 1)2, A, 4  " 

C        FIRST  TRIANGULAR  TEST  COMPLETED  

'4   TRI  2  ■'=  D  +  C-E 

IF(TRI2)2,5,5 
5" "TRI  2  =  ~"D-C'*E 

IF(TRI2)2,6,6 
6   TRI 2  =  -D+C+E 
IF(TRI2)2,7,7 


SECOND  TRI  ANGULAR~7EST~  COMPLETED" 

7  TRI3  =  A+C^F 

IF (TRI 3) 2, 8, 8 

8  TRI3=A-C+F 

IF(TRI3)2,9,9 

9 TRI  3  =  -A  +  C+F 

~  IF(TRI3)2,  10, 10 

THIRD  TRIANGULAR  TEST  COMPLETED 
i'O  "TRI  4  =  D  +  B-F 

IF  (TRI  4)  2",  11,11 
IT  TRI  4  =  D-B>F 

IF(TRI4)2,12,12  

"17   TRI4  =  -D  +  B  +  F 

IF(TRI4)2, 13, 13 
FOURTH  TRIANGULAR  TEST  COMPLETED 
13   DELI  =  DELTAIA.B.E) 
"  DEL2"=  OELTA(D,C,E) 

DfL3  =  DELTA (D,B,F) 

DEL4  =  DELTAl  A,C,F) 
DELX  =  DEL1*DEL2*DEL3*DEL4 
~"N  =  A  +  B  +  C  +  D 

PHZ  =  (-1.0)*»N 
"SUM  =  0.0 

_AK  =  0.0      

"24   SI"  =  A+H-E-AK- 

IF(S1)22,16, 16 
16 "  M  =  SI 

FS1  =  FACT(M) 
S2  =  C+D-E-AK 
IF(S2)22,i8, 18 
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18  M  =  S2 

FS2  =  FACT(M) 

53  =  A  +  OF-AK  

IF (S3) 22, 19, 19" 

19  M  =  S3 

FS3  =  FACT(M) 

54  =  D+B-F-AK 
IF(S4)22,20,20 

20  M  =  S4 

FS4  =  FACT(M) 

55  =  -A-D+E+F+AK 
IF(S3) 17,21,21 

21  M  =  S5 

FS5  =  FACT(M) 

56  =  -B-OE  +  F  +  AK  

TF(S6) 17,2  3,23 

23   M  =  S6 

FS6  =  FACT(M) 
N  =  AK 
"'SPHZ  ="  (-1.0)  "N" 

TOP  =  A»B*C+U+1,0-AK 

~M  V'TOP 

FTOP  =  FACT(M) 
"TOP  =  SPHZ»FTOP 
M  =  AK 
FAK  =  FACT(M) 

DENOM_=  FAK»FS1*FS2«FS3«FS4»FS5»FS6 

SUM  =  SUM  ♦  (  TOP/OENOM)" 
17   AK  =  AK  ♦  1.0 
"GO  TO  24 

22  SIXJ  =  PHZ*DELX»SUM 
"RETURN" 

END 
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DELTA.  .. 

FUNCTION  DELTA(A,B,C) 
ROTENBERG  ET.  AL.  PAGE  13  (2.4) 

FUNCTION  DELTA  (  A'/BVCT  " 

SI  =  A  +  B  -  C 

M  =  SI 

FS1  =  FACT(M) 
"  S2  =  A  ♦  C  -  B 

M  =  S2     ' 

~FS2  "=  FACT(M) 

S3  =  B  +  C  -  A 
"~M  =  S3 

FS3  =  FACT(M) 
""OENOM  =  A+B+C+1.0 

M    =    DENOM 
T0~~=    FACT(M) 

DELTA    =       SQRT( (FS1»FS2*FS3)/FD) 
"RETURN" 

END 
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F ACT • • • 

CALCULATES    N    FACTORIAL 

FUNCTION    FACT(N) 

A  =  N 

FACT=0.0 

IF(A.LT.O.O)     RETURN 

FACT=1.0 

IF(A.EO.O.O)     RETURN 

FACT=FACT*A 

00  T~T"»V»100~ 

B=I 

C  =  A-8 

IF(C.EQ.O.O)     RETURN 

FACT=FACT»C " 

CONTINUE  

MR  I T  ET671 00 ) 


100 


FORMAT! 19H0FACT0RIAL 

STOP" 

END 


OVERFLOW 


-103- 
QTAVE 

The  required  card  decks  are: 


1. 

SEARCH 

2. 

COREL 

3. 

OUTPUT 

4. 

DERIV  3 

5. 

JPLOT 

6. 

CEFMAT 

7. 

THEORY  6 

8. 

QSQUAR 

9. 

HERM 

see  listings  for  OPTIC 
see  SHARE  no.    884 


The  listings  of  decks   1-7  follow, 


-  1  04- 

SEARCH.. . 

THIS  PROGRAM  USES  THE  METHOD  OF  DIRECT  SEARCH  TO  OBTAIN  BEST 

VALUE  OF  CHI»«2 

INPUT  DATA 

NE=NO.  OF  ENERGY  LEVELS  MEASURED 
NP=NO.  OF  SHIELDING  PARAMETERS 

NPP=NO.  OF  SHIELOING  PARAMETERS  WITH  NON-ZERO  STEP  SIZE 
~ITIME=  MAX.  XEQ  TIME  ALLOWED 

E( I  )=MEASURED  ENERGY  LEVELS  

~R(I)=ERROR  IN  MEASURED  ENERGY  LEVELS" 
V(I)=LEVEL  IDENTIFICATION 
P(2,  I  )  =  INITIAL  VALUE  OF  PARAMETERS 

X(I)=STEP  SIZE  FOR  EACH  PARAMETER  FOR  INITIAL  CYCLE 
XMIN=2**N  ,  WHERE  N  IS  MAX  NO.  OF  CYCLES  ALLOWED  I.E.  FIMAL 
STEP  SIZE  WILL  BE  X( I  )/2*»N 
DIMENSION  E(  100)  ,R(100),V(  100)  ,  P  (  3V20T."X  (  20  JVQSQ  (  3  )  ,  DP  (20  )  .SC1500) 
DIMENSION  AB(20) 
DIMENSION  GIVE( 300)  * 

COMMON  GIVE,S,E,R,V,P,X,QSQ,NE_,NP,  IQ , JQ,KQ, I  OP, XTEST ,NFORM , NMGVE , A 
""COMMON  /COMTIM/ITIME 

READ   (5,99)NE  ,NP,NPJS  ITJME 

99  FURMAT( I3/I3/I3/I6) 
CALL  ICLOCKJ INTIME) 
~I  TIMEMNTIME  +  ITIME 

READ   (5,800) (E(  I  )  ,R(  I  ) ,V(  I ) ,  1  =  1, NE) 
800  FORMAT ( 3F20.5)' 

READ   (5,50  0)AA,BB, ( A8 ( I  ) ,  1  =  1, NP) 

"500  FORMAT  (  12A6) 
NFORM=0 


1000  READ   (5,200)  (  P  <  2,  I)  ,X(  I  ),"  I  =  1,NPT 
200  FORMAT(2F20.5) 

"READ   (5,300)XMIN 

_300_F0RMAT  (  F_20  .b) 

NMOVE=0" 

0S0(2)=QSQUAR( 1,2,2) 

NF0RM=1 

A=NE-NPP 

CHI=QSQ(2)/A 

WRITE  (6,700)  (AB(  I  )  ,  1  =  1, NP) 

700"  FORMAT  (  1  H  l,  (  1  OX  ,  A6  ,  1  4X  ,  A6,  14X  ,  A6  ,  1  <*X  ,  A6  ,  1  4X  ,  A6  ,  1  4X  ,  A6  )  ) 

WRITE  (6,701) 

701  F0RMAT(9X,7H  CHI»«2//29H  INITIAL  VALUES  OF  PARAMETERS) 
WRITE  (6,702)  (P(2, I  ),  1  =  1, NP) 

""702  FORMAT  (6E20.8) 

WRITE  (61703)CHI ,  (X(  I  )  ,  I=1,NP) 

"703  FORMAT  (  1PE20.7//2AH  IN  ITT  AL'"VALUE"S"  OF  STEPS/ f  0P6E20.  8  )  ) 
C  FORMING  PATTERN 

XTEST=1.0 
13  NZERO=0 

NF0RM=NF0RM+1 
DO  2  1=1, NP 


105- 


50 


"2  0 
11 

10 
~\2 


100 
101 


8 
18 

14 


15 
16 

17 


IF(X(  I  )  ) 

CONTINUE 
P( 1, I )=P 
QSQ(  1) =0 
P( 3, I )=P 
0SQ(3)=Q 
IF (0SQ(2 
IF(QS0(2 
DPI  I )=-l 

"QSQ(2)=Q 
GO  TO  7 
DPI  I  )=0. 
NZERO=NZ 
GO  TO  7 
DPI  I  )  =  +  l 

'0SQ("2)=0 
Pi2f  I )=P 
CONTINUE 
CHEC 
IFIN2ER0 
XTEST=2. 

"CALL  I  CL 
IF (LPTIM 
IF  (XTEST 
CALL  OUT 
GO  TO  10 
DO  12  1= 
XI  I  )=X'(  I" 
CHI=OSQ( 
WRITE 
FORMAT (6 
WRITE 
FORMAT! 1 
GO  TO  13 
MAKE 
NMOVE=0 
NMOVE=NM 
DO  14  1= 
PI  I,  I )=P 
QSQI  1)=Q 
IF (QS0I2 
CSQ(2)=Q 
DO  16  1= 
PI2, I)=P 
CALL  I CL 
IFILPT  IM 
GO  TO  18 
END 


50,5,50 

(2,  I  )-X(  I )    _ 
SOUARC 1,1,2) 
(2, I  )+X(  I  ) 
SQUARI 1,3,2) 
)-OSQ« 1) )3,3,4 
)-QSQ(3) )5,5,6 
.0 
SQU) 


0 
ERO+1 

.  0 

SQ  (  3  ) 

(2,1  )+DPI I  )»X(  I  ) 

K  TO  SEE  IF  PATTERN  IS  NON-ZERO 
-NP)8,9,8 

0»XTEST     

OCKILPTIME) 

E.GT.ITIME)  GO  TO  11 

-XMIN) 10, 10, 11 

PUT 

00 

l.NP 

)/2.0 
2)/A 

(6, 100) (PI2.I ) , 1=1, NP) 
E20.8) 

(6,101)CHI,NFORM,NMOV£,XTEST 
PE20.7/15X,7H  NFORM=  I6,7X,7H  NMOVE=  16, 7H  XTEST=  0PF9.1/) 


PATTERN  MOVES 


OVE-t- 
1,NP 
(2,1 
SQUA 
)-QS 
SOI  1 
1,NP 
11,1 
OCKI 
E.GT 


1 


)+DP( I )»Xl I ) 

R(  1,1.1") 

Q( 1 )  )13,  13,15 

) 


) 

LPTIME) 

.  ITIMETXO  TO  11 


-106. 


TO 


11 


200 


12 


100 
90 


COREL. 
THIS  S 
SQUARE 

"SUBROU 
DIMENS 
DIMENS 
DIMENS 

" COMMON 

COMMON 

~NK  =  NE~ 

DO  5  I 

"DO  5  K 
NE=K 
CALL  I 
IF(LPT 
5~DER( I, 
DO  10 

"DO  10 

AINV( I 
"DO  10 
AINV ( I 
CONTIN 
NPM1=N 
DO  11 
IP1=I+ 
DO  11 
AINV( I 
CONTIN 
WRITE! 
FORMAT 
CALL  M 
DO  12 
DO  12 
A  (  I  ,  J  ) 
CONTIN 
WRITE 
FORMAT 
RETURN 
END 


URROUTINE  CALCULATES  THE  CORRELATION  MATRIX  FOR  A  LEAST 

ANALYSI S. 
TINE  COREL 

ION  A( 5,5) ,DER(5,20), AINVI  50, 5 ) , B { 50, 1 ) 

ION  E( 100) ,R( 100), V( 100) ,P{3,20) , XI20) ,0SQ(3) , DP ( 20 ) , S ( I  50  0 ) 
ION  GIVEOOO) 

GIVE,S,E,R,V,P,X,QSQ,NE,NP,IQ, JQ, KQ , I  OP ,XTEST , NFORM , NMOVE 

/COMTIM/ITIME 


=  1,NP 
=  1,NK 


CLOC 
IME. 
K)=D 
1  =  1, 
J=i, 
,J)  = 

K=r, 
,jj  = 

UE~ 

P-l 

1  =  1, 

1 

J=IP 

,J)  = 

UE~ 

6,20 

(36H 

ATIN 

1  =  1, 

J=l, 

=  SQR 

UE 

(6,  1 

(44H 


GO  TO  90 


K(LPTIME) 
GT. ITIME  ) 
ERIV(TVK)— 

NP 

0,0 

AINV(  I  , J)+DER(  I,K)»DER(J,K)/(R(K)*R(K)) 


NPM1 


1,NP 
AINVU,  I  ) 


0)  (  (  AINVU  ,  J  ),  1  =  1, NP),  J=1,NP) 

INVERSE  "OF  BEST  CORRELATION  MATR I X// ( 5E20 . 8 )  ) 

V( AINV,NP,B,0,DETERM) 

NP 

NP 

T(ATNV'(T,J)  ) 


00) ( ( A( I , J),  1=1, NP) ,  J=1,NP) 

BEST  CORRELATION  MATRIX  (SORT  OF 


ELEMENTS)//(5E20.8) ) 


_J  07  _      , 

C      OUTPUT... 
"C      THIS  SUBROUTINE  CONTROLS  THE  OUTPUT  PHASE  OF  SEARCH 

SUBROUTINE  OUTPUT 

D I  MENS  I  UN  E(IOO) ,R( 100),V( iOO) , P ( 3, 20 ) , X ( 20 ) ,QSQ ( 3 ) , DP { 20 ) , S ( 1500 ) 
DIMENSIUN  GIVE(300) 

COMMON  GIVE,S,E,R,V,P,X,QSQ,NE,NP, 10 , JQ.KQ, I  OP ,  XT EST , NFORM,NMQVE , A 
CHI=QSQ(2)/A 
"XTEST=XTEST/2.0 

WRITE  (6,200)  

~200~  FORMAT  (38H  FINAL~VALUES  OF  PARAMETERS~AND~CHl  »«2 T~ 
WRITE  (6,800) (P(2, I  >  ,I  =  1,NP) 

800  FORMAT(6E20.8) 

WRITE  (6, 801 ) CHI , N FORM , NMOVE , X TEST 

801  FORMAT{ 1PE20.7/15X.7H  NFORM=  I6,7X,7H  NMOVE=  16, 7H  XTEST=  0PF9.1/) 
WRITE  (6,300) (X(  I  )  ,  1  =  1, NP) 

300"  FORMAT  (  2 2HFlNAL"VALUFS0F~STEPSYr6F20T8TJ 

CALL  CCKEL 

CALL  THEORYfl,  2,2,0) 

RETURN 

ENO 


-. 


-  108 


TO" 

20 
30 


D  E 
THIS 
OF  T 
POIN 
DERI 
ONLY 
FUNC 
"COMM 
IF  (  I 
DERI 
GO  T 
DERI 
RETU 
END" 


R  I  V 

FUNCT 
HE  SUB 
T  JP. 
V3  IS 

IN  TH 
TION  0 
ON/COM 
P-1J10 
V=DR(J 
0  30 
V=OSU 
RN 


ION  CALCULATES  THE  PARTIAL  DERIVATIVE  OF  THE  OUTPUT 
ROUTINE  THEORY  WITH  RESPECT  TO  THE  PARAMETER  IP  AT  THE 


FOR  THE  SPECIAL  CASE  OF 
E  CASE  OF  C2  SYMMETRY. 
ERIVI IP.JP) 
DER/DR(300),DS(300) 

,10,20 

PI 


THE  SHIELDING  PARAMETERS 


P) 


109 


JPLO 

ON  T 
AND 


1000 

999 

1001 


980 


994 
A3 


44 
45 


46 


952 
991 
997 

950 

995 


JPLOT. 
T,  THI 
HE  OFF 
THE  AB 
S  =  ARR 
NP=  NO 
NL=  NO 
X0=  LO 
Xl=  UP 
X=~UNI 
DD=  AN 
LAB  =  N 
WHEN  L 
CONSIS 
PLUS  S 
SUBROU 
DIMENS 
COMMON 
IF(LAB 
READ 
FORMAT 
JD=NP/ 
X2=(X1 
WRI  TE 
FORMAT 
15X,  10 
I  I  I»I  I 
DO  998 
ABC=X* 
DO  9  94 
AB(  I  )  = 
IF(S(J 
AB(  I  )  = 
AB(2)= 
GO  TO 
IF (S( J 
AB( 16) 
AB(  17) 
GO  TO 
PS=100 
IPS=PS 
ITEST= 
DO  997 
IF{ IPS 
IF( IPS 
ITEST= 
CONTIN 
IDEL=I 
GO  TO 
AB(  I  )  = 
GO  TO 


S  SUBROU 
-LINE  PR 
SCISSAS 
AY  TO  BE 
.  OF  ELE 
.  OF  LIN 
WER  LIMI 
PtzR  LIMI 
TS  PER  L 

ALPHAME 
0.  OF  TH 
AB=l  THE 
T  OF  7  W 
IGN. 

TINE  JPL 
ION  S(15 

GIVE 
-1)  1000, 
(5,999) ( 
(9A6) 
NL 

+X0)/2.0 
(6,980)0 
(18H1  RE 
2\\  »I  I  I  I 
IIIIIIT* 

J=1,NP, 
FLOAT( J) 

I  =  1 ,  1 7 
AA(7) 
)-X0)43, 
AA(8) 


TINE  PLOTS  A  SINGLE  VALUED  FUNCTION 

INTER.   THE  ORDINATES  ARE  PLOTTED  ACROSS  THE  PAGE 

ARE  PLUTTED  DOWN  THE  PAGE 

PLOTTED  (ONE  DIMENSIONAL) 
MENTS  IN  S 

ES  THAT  PLOT  IS  TO  COVER 
T  OF  ORDINATES 

T  OF  ORDINATES 

INE  (SCALE  FACTOR) 

RIC  TITLE  FOR  PLOT,  10A6 

E  PLOT,  MUST  START  WITH  1 

PLOT  SYMBOLS  ARE  READ.   THE  PLOT  SYMBOLS,  AA( I ) , 
ORDS"  THAT  CONTAIN  6  BLANKS  OR  5  BLANKS  AND  ONt 


OT(S,NP, 
00)  ,GIVE 

,  S 
1000, 100 
AA(I),  1=1,9) 


NL",  XO  ,"X  1  ,  X',  X  I  ,  00V  LAB  T 
(300) ,AA(9) ,AB(20) ,DD( 10) 


L 


D,X0,X2, 
SULTS  OF 
I  I  I  I  I»I  I 
I  I  II  I  III 
JO 
♦  XI 


XI 

JPL0T//10A6//12X,F7.2,43X,F7.2,43X,F7.2/ 
I  I  I  I  I  I  I  *  I  I  I  I  1  I  I  I  I  *  I  I  I  I  I  I  1  I  I  *T  I  I  I  1 1  I  I  I  •  1 1  I  I  1  I 
I »  1 1  I  I  I  I  I  I  1*1  I  I  I  I  I  I  I  !•  ) 


44,44 


AA(9) 

990 

)-Xl  )46, 46,45" 

=AA(8) 

=AA(9) 

990 

. 0  * ( S ( J ) -xo ) / ( X 1 -X 0 ) +075 


0 

1  =  1 

-ITL 

-6-1 

ITES 

UE 

PS-I 

(1,2 

AA{  1 

990 


,17 

ST  )952,995,9 

TEST)950,99l 

T  +  6 


52 

,991 


TEST 

,3,4,5), IDEL 
) 


110- 


4  AB( I )=AA( b) 
GO  TO  990 

3  AB(I)=AA(4) 

GO  TO  99  0 
2  AB(  I )=AA(3) 
' GO  TO  990 
1  A8(  I  )  =  AA(2) 

GO  TO  990 

5  AB( I  )  =  AA(6) 

"990  VJRITE  T6~,  989TABC 7TABTKT,  K=l,  17) 
989  FORMAT (7X,F6. 1,3H   I    17A6) 
998  CONTINUE 

RETURN 

END 


-Ill 


C  E  F  M  A  T 

THIS  SU6R0UT 

EARTH  IGN  IN 

HIGHER  J  MAY 

ARE  CHANGED. 

INPUT  D 

P(2,I)= 

C4-2,C4 

NTERM=N 

"NJ  =  NO.~ 

XJ=MIN. 

I  W  =  0'UMM 

IDEG=l* 

THET4(N 

IDENTIF 

Ajn)=j 

NOP=NO. 
~  EACH  J 
OA,  OB, 
KTH'J  V 

GO  FROM 

FOR'  MJ  = 

SUBROUTINE  C 

"DIMENSION  OA 

DIMENSION  P( 

COMMON/COMTH 

READ(5,699) ( 

"699  F0RMAT(6E12. 

1  READ   (5,100 

"  100  FORMAT(2I2,F 

DU  999  M=1,N 

READ   (5,101 

101  F0RMATOE20. 
~999  CONTINUE 

DO  1000  1=1, 
DO  1000  J=l, 
DO  1000  K=l, 
0A( I  ,J,K)=0. 
_0B( I  ,  J,K)=0. 
OC  (  I  ,J,K)=0. 

1000  CONTINUE 

DO  998  K=1,N 
READ   (5,103 

"103  FORMAT  (12) 
DO  997  L=1,N 
READ   (5,  102 

102  FORMAT  (212, 
997  CONTINUE 

AK  =  K 

LMAX=AK+XJ+0.52 
DO  996  L=1,LMAX 


INE  C 

NTER 

BE  H 

THE 

ATA 

CEF  P 

4,C4- 

0.  OF 

OF  DI 

J  VA 

V  VAR 

FOR  I 

,M)  =  R 

ICATI 

VALU 

OF  C 

VALUE 

OC(K 

ALUE 

1  TO 
J 

EFMAT 
(8,  16 
3,20) 
/CEF( 
P(2,  I 
5) 

)NTER 
10.5, 
TERM 
)  (THE 
5,F10 


ALCUL 

M  TER 

ANDLE 

ELEM 

ARAME 
4,C60 

OPTI 
FFERE 
LUE  - 
IABLE 
NTEGR 
EDUCE 
ON 

E  FOR 
ARDS 

DIVI 
,I,J) 
FOR  A 

( 10,  J 
,16), 

13,  13 
),  1  = 


ATES  THE  CEF  MATRIX  ELEMENTS  OF  A  RARE 
MS  OF  THE  OPTICAL  SPECTRA.   THE  MAX  J=lt»/2. 
D  IF  THE  APPROPRIATE  U  I  MENS  I  ON  "STATE  MtN  T S 
ENTS  ARE  CALCULATED  FROM  GIVEN  CEF  PARAMET. 


TERS 
tC62, 
CAL  L 

NT  J  VALUES 
1.0 


IN  THE  FOLLOWING  ORDER  C20 , C22 , C40, C42 , 

C6-2,C64,C6-4,C66,C6-6 

SJ  TERMS 

IN  TERMS^   " 


AL  J, 

D  MAT 


THE 
WITH 
DED  B 
=  IJ  0 
LPHA, 

•  \  y 

0,KQ) 
0B(8, 

»2) 
1,  14) 


2  FOR  HALF  INTEGRAL  J 
RIX  ELEMENTS,  N  I"S  ORDER, M  IS  TERM 

ITH  TERM 

OPERATOR  EQUIVALENT  MATRIX  ELEMENTS  FOR 

Y  3 

PERATOR  EQUIVALENT  MATRIX  ELEMENT  FOR 

BETA,  GAMMA  RESPECTIVELY,  '  I  AND  J 
S  INDEX  FOR  MJ=-J  AND  2J+1  IS  INDEX 


16,  16), OC (8, 16, 16) ,THETA(3,20) 
ANM(3,4,2) ,AJ(20) ,SIGN(2) 


M,NJ,XJ,I W, IDEG 
212) 


TA(N,M) ,N=1,3),AJ(M) 
.5) 


8 

16 

16 


J 
)NOP 

OP 

) I , J, (OA(K 

3F20.5) 


I  ,J),OB(K, I ,J) ,OC(K,I,J)  ) 


112 


996 


995 
"998 


"1001 
2 


. 


11  1 

99  3 

99  2" 

991 
994 


903 


MIN=l+ 
JMAX=2 
M  A  X  =  J  M 

DO  996 
NA  =  JMA 
OA(K,M 
OBIK.M 

OC  |K,M 
CONTIN 
DO  995 
DO  995 
OA (K,N 

0  H  (  K  ,  N 
0C(K,N 
CONTIN 
CONTIN 

DO  100 
DO  100 
DO  100 
ANM( I , 
CONTIN 
IA  =  0 

is=r 

DO  994 
MAX=N+ 
DO  994 
DO  994 
IFIL-N 
CONT IN 
IF (M-L 

1  A  = I A  + 
IF( I A- 
ANM(N," 
GO  TO 
ANM(N, 
CONT IN 

SIGN(  1 
"SIGN  (2 
1  =  1 

I  J  =  AJ( 
JMAX=2 
DO  903 
DO  903 
DO  903 
CEF( J, 
CONTIN 
DO  901 
DO  901 
IF( J-K 


1 

.0* 
AX  + 
N  = 
X+l 
AX, 
AX, 
AX, 
OE 
L  = 
N  = 
»L) 
,L) 
,L) 
OE 
OE 

1  I 
1  J 
1  K 
J,K 
UE~~ 


( AK+XJ ) ♦ 1 . 01 
1-L 

MIN,MAX 

-N 

N)=OA(K,NA,L) 

N)=OB(K,NA,L) 

N)=OC(K,NA, L) 


1,  JMAX  " 

1,L 

=OA(K,L,N) 

=08(K,L,N) 

=OC(K,L,N) 


SETTING 
=  U3 
=  1,4 
=  1,2 
)=0. 


OP  OF  ANM 


1 


N-1,3 

M=  1  ,  M  A_X_ 

L  =  l  ,Z 
)  111,  111,994 
OE 

)994,993,993 
1 

10)992,991,992 

M,L)"=P(KQ,"IA) 

994 

K,L)=P( JO, IA) 

OE 

CALCULATE  CEF  MATRIX  ELEMENTS 
)  =  1  .  0 

)=-~iYo 


I )-XJ+.Oi 
•0»AJ( I )+1.01 

J=1,JMAX 

K=l , JMAX 

L=i,2 
K,L  )=0.0 
UE 

J=1,JMAX 

K=1,J 
)800,801,800 


113 


801 


80 

"80 

"  8C 
80 

80 
80 
81 


"89 
90 


1002 


M=i 

GO  TO 
IF< J-K 
M=2 
GO  TO 
IF (J-K 
M  =  3 
GO  TO 
IF(  J-K 

CONTIN 

DO  899 

CEF( J, 

1     08 

CEF(Kf 

"CONTIN 

CONTIN 

WRITE( 

FORMAT 

1R2=F10 

1  FiO 

2T377H 

38H   B6 

47H  8  66 

RETURN 

"END 


810 
-2)901,802,803 

810 
-4)901,804,305 

810 
-6)901,806,901 


UE 

L  = 
K,L 
(IJ 
J,L 
UE 
UE 
6,1 
(32 
.3, 


1,2 

)=SIGN(L)«(THETA(l,I)«0A(lJ,J,K)»ANM(l,MfL)+THETA(2,I)» 
, J,K)«ANM(2,M,L )+THETA( 3, I)*OC( I  J, J,K)»ANM( 3,M,L )  ) 
)=SIGN(L)»CEF(J,K,L) 


002)  (P(2,I  ),  1  =  1,  15) 

HI   VALUES  OF  CEF  PARAMETERS  USED/  8H   A20R2=F 10. 3 , 7H  A22 

7H  B22R2=F10.3,7H  A40R4=F10 . 3 , 7H  A42R4=F 1 0. 3/8H   B't2R4  = 


A44R4~=F10.3,7H  B44R4=F  10.  3  ,7H  ~A60R6  =  F10.3  ,7H_A62R6  =  F13.  3  A 
2R6=F10.3,7H  A64R6=F 10. 3 , 7H  B64R6=F 1 0. 3 , 7H  A66R6=F10.3, 
R6=F10.3) 


14 


THE 

THIS  S 

__  CALCUL 

~1N  TER 

I 

N 

I 

D 

P 

S 

C 

D 

Y 

" SURROU 

COMMON 

COMMON 

COMMON 

"COMMON 

D I  MENS 

DIMENS 

DIMENS 

"DIMENS 

COMPLE 

"COMPLE 

IFINFO 

2CALL  C 

DO  3  I 

~~WTI  )=C 

3  CONTIN 
1  =  1 

DO  4  J 
JP1=J+ 
DO  4  K 
H( I )=C 
1  =  1  +  1 

4  CONTIN 
SPARE( 
DELTA= 

CALL  H 

"MULTIP 
DO  100 
TEMP=R 
IF (TEM 
CTEMP= 
GO  TO 

1001~CTEMP= 

1002  00  100 

1003  C ( I , J) 
1000  CONTIN 

REARRA 
DO  5  I 


0  R  Y  6  . 
UBROUTINE 
ATEStDELT 

MS  OF  SHI 
NPUT  DATA 
P0INT=N0. 
F  NPOINT 
EPENDENT 
LOTTER 
CALt=SCAL 
20,C22=CE 
D=LAhEL  F 
MIN, YMAX= 
TINE  THEO 
GIVgiS.E 
/COMDER/D 
/COMHER/E 
/COMTH/CE 
ION  QZZ(  1 
ION  GIVE( 
ION  E(  100 
ION  EVECI 
X  H,C,CTE 
X  CSQRT, 
RM)  1,2,  I 
EFMATI  1,2 
=  1,13 
EF(  I,  I,  1) 
UE 


DIAGONAL 
A  E)T,  TH 
ELDING  PA 

OF  POINT 
IS  GREATE 
QUAD  SPLI 

E  FACTOR 
F  PARAMET 
OR  PLOTS 
MIN.  AND 
RY( 10, JQ, 
,R,V,P,X, 
R(300),DS 
VENUP,W(  1 
F( 13,  13,2 
3) ,QXY(  13 
300) 

)  ,R( 100), 
(  13) ,EVEC 
MP, CMPLX 
CONJG,  CA 


IZES  CEF  MATRIX  PROVIDED  BY  CEFMAT  AND 
E  TEMPERATURE  DEPENDENT  QUAD  SPLITTING 
RAMETERS. 


S  THAT  ARE  TO  BE  CALCULATED  FOR  PLOTTING 
R  THAN  OR  EQUAL  TO  150  THE  TEMPERATURE 
TTING  IS  ALSO  PLOTTED  ON  THE  MOSELY 


FOR  TEMPERATURE  SCALE 
ERS 

MAX.  VALUES  OF  QUAD  SCALE  IN  CM/SEC 

KQ, IOP) 

QSQ,NE,NP, IQ,JQ,KQ,IOP,XTEST,NFORM,NMOVE,A 

(300) 

3)  ,H(78)  ,SPARE(79) 

) 

It  00 (10) ,C( 13, 13), EX (13, 100) 

V(100),P(3,20),X(20)  ,frSQ(3),DP(20)  ,S(15Q0) 
R  (  1  3  ) 


BS 


,2) 


=  1,12 

1 

=JP1,13 

MPL  X ( C  EF  (  J  ,  KVHTCE  F  (  J  ,  K  ,  2  )  ) 


UE 

79)=0.0 

1.0E-10 

ERM(H,  13, CO, DELTA,  IT) 

LY  C  BY  PROPER  PHASE  FACTOR 

0  1=1,13 

EAL(C( 1,7)) 

P.NE.O.O)  GO  TO  1001 

CSQRTI-CONJGICU  ,12)  )/C(  1,2)  ) 

1002 

CONJG(C( 1,7) )/CABS(C(I,7))_ 

3  J=i, 13 

=CTEMP«C( I ,J) 

UE 

NGE  ENERGIES  AND  WAVE  FUNCTIONS 

=  1,12 


115- 


IP1=I+1 

DO  5  J  =  I PI  ,13 
_IF(W(I)-W(J))5  t2.lL 
1    TEMP  =  W( I  ) 

W  (  I  )  =  W  (  J  ) 

W( J) =  TEMP 

DO  8  K  =  l,  13 

CT EMP=C( I ,K) 

C( I,K)=C( J,K)  

'8  C(  J,K)=CTEMP~ 
5  CONTINUE 


RENORMALIZE  THE  ENERGIES 

TEMP  =  W(  1  ) 
"00  9  1=1,13 
9  Wl  I  )=W( I  )-TEMP 
"CALCULATE  MATRIX~£LEMENTS"OF  "3JZ»«2-jrj+n~AN0~3A2nP»»2  +  JM«*21 

DO  10  1=1,13 
~QZZ(  I )=0.0  " 

DO  10  J=l,13 
"OZZ(  I  )=QZZ  (  I  )*REAL(C'ONJG(C'(I',"J))»"C(r,J))«(3"."0»(FL'OAT{J)-r.O)" 

1  »«2-42.0) 

l'O  CONTINUE 

DO  1  1  I  "=  1,13 
"OXY( I )=0.0 

DO  11  J=3,13 
" AJ=J-7 

_OXY(  I  )=QXY(  I  )*1.5»REAL(C0NJG(C(  I  ,  J  )  )  »C  (  I  , J-2)*C(  I  ,  J  )j»_ 

1  CONJG(C(  I  ,J-2  )  )  )*SQRT(  (  6."0'*A J ")  *  (  5.  0'*A  J  )  •  (  7.  0- AJ  )  •  (  8.  0-A  J  )  ) 
11  CONTINUE 

RE  AD  (5,  100)  NPOINT,  SCALE  ,C  20,  C  22 
100  FORMAT( I4,F10.5/2F20.5) 
"BSD.  695056 

E20=0. 769854  

~~HCE  =  0.8610  308£-~5 

R2=0.19 
" R3=75.5 

ALPHA=l.0196651E-02 

CALCULATE  "THE  PARTITION  FUNCTION 

1_  I  F  (  I  OP  )  1 2  »_13'JJ 

13  M1=NP0INT 

GO  TO  15 
i"2"Ml=NE 

IF(NFORK)  15,  U,  15 
"DO  19  M=1,M1 

T  -  V  (  M  )  _____ 

00  19  N=l, 13 

EX(N,M)=EXP(-W(N)/(B*T)  ) 

CONTINUE 

DIR1=0.5*E2Q*4.0«HCE»C20/R2«(-1.0) 

DIR2=0.5*E2Q»4.0»HCE*C22/R2»l-1.0) 

C4F=0.5*E2Q»ALPHA»R3»(-1.0) 


- 

LA 

16 

18 

19 

116 


15 


50 

51 
52 


71 


12 


7  0 


21 


22 


000 


001 


002 
2  3 


20 

25 

200 


CONT 
FEED 
GAMM 
GAMM 
GO  T 
GAMM 
GO  T 
GAMM 
CONT 
DO  2 
IF(  I 
T  =  SC 
DO  7 
EX(  I 
JAP  = 
I  F  (  I 
PART 
DO  2 
PART 
OZZA 
QXYA 
DO  2 
QZZA 
QXYA 
SA  =  C 
SAA  = 
IF(D 
S(  J) 
GO  T 
SB  =  C 
SBB  = 
S(  J) 
IF  (  I 
CONT 
DR(  J 
DS(  J 
CONT 
IF(  I 
WRIT 
FORM 


31 

201 
30 


1 

2117H 

3 
DO  3 
DO  3 

'  EVEC 
EVEC 
WR  IT 
FORM 
CONT 
WRIT 


INUE 

IN 
AR  =  P 
AS  =  P 
0  (5 
AR  =  P 
0  52 
AS  =  P 
INUE 
0  J  = 
OP)  7 
ALE* 
2  1  = 
,1)  = 
1 

OP.N 
=  0.0 

1  1  = 

=  PAR 

=  0.0 
=  0.0 
2  1  = 

=  ozz 

=  OXY 
4F*Q 
SA*( 
IR2) 
=  SAA 
0  20 
4F*Q 
SB»( 
=  SQR 
0P)2 
INUE 
)=-( 
)=(S 
INUE 
0P)2 
E(6, 
AT(3 
4 
( 
0 

0  1  = 

1  J  = 


THE  PARAMETERS 

(KQ,1)    

(KQ,2) 
0,51 ) , IQ 
( JQ, 1) 


AND  CALCULATE  THE  QUAD  SPLITTING 


( JQ,2) 


1,M1 

0,71,70 

FtOATI J) 

1,  13 

EXP(-W(  I )/(B*T)) 


ETOF"  JAP  =  J 

1,13 

T*EX(  I  ,  JAP) 


1,13 

A  +  OZZ(  I  )*EX(  I, JAP) 
A+QXY(  I  )*EX(  I, JAP) 
ZZA/PART 

1.0-GAMMAR)+DIRl»GAMMAS 
2001,2000,2001 


02 

XYA/PART 

1.0-GAMMAR)+DIR2*GAMMAS 

T(SAA**2+1.0/3.0*SBB»»2) 

3,20,23 


3AA»SA+SBB«SB/3.0)/S 
AA»DIRl+SBB»DIR2/3.0)/S 


4,25,24 

200)  IT 

2H       NUMBER" 'OF  ITERATIONS 

8H0   ENERGY 

CM-1)    MJ=     -6 

+  1  *2 
1,13 
1,13 
=  REA 
=  AIM 
201) 
10.3 


-1 


R(  J) 
I  (J) 
E{6, 
AT(F 
INUE 
E(6,202)QZZ,OXY 


IN  HERM     15   / 

MIXING  COEFFICIENTS/ 
-5       -4       -3       -2 
+  3       +  4       *5       +  6// 


L(C(  I, J)  ) 

AG(C( I ,J) ) 

W(  I  ), (EVECRt J) ,  J=l,13),  (EVEC I (K) , 

,6X  13F8.3/16X  13F8.3) 


K=l,13) 


1  17 


2C2  FORM 

1  ♦  J 

READ 

"203  FORM 

WRI  I 

204  FORM 

1  FIO 

I  F  (  N 

CALL 

CALL 

RETU 

END 


32 
24 


AT(18H0   ( 3JZ*»2-J(J+1 ) )//6E20.3/6E20.8/E20.8/22H0   l.5»(J*»*2 
-••2)//(6E20.8) ) 

(5,203)  DD,YMIN,YMAX 

AT  (  10A6/2F20.5) 

E(6,204)  DD, SCALE,  (S( I  )  ,  1  =  1, Mi) 

AT(22H1   RESULTS  OF  PLOT  OF  10A6/10X  15H  SCALE  FACT3R= 

.2//I2X  20F6.3) ) 

POINT. LT. 150)  GO  TO  32 

CPLOT ( S,M1 ,M1 ,DD, YMAX,YMIN,  1  ) 

J  PLOT  (S,  Ml,  120,YMIN,YMAX,  SCALE",  0.0 ,DD ,  1  T 
RN 


-1  18- 
HITEMP 
The  required  card  decks  are: 


1. 

THEORY  5 

2. 

DERIV    2 

3. 

SEARCH 

4. 

QSQUAR 

5. 

OUTPUT 

6. 

COREL 

7. 

JPLOT 

see  listings  for  QTAVE 


The  listings  of  decks    1-2  follow, 


-119- 

C         THEORY   5... 

C      THIS  SUBROUTINE  CALCULATES  THE  QUADRUPOLE  SPLITTING  IN  THE  FORM 

C _0_F_A/T*B  

C   '  I NPUT  DATA  ' 

C  TEMP='INITIAL  TEMPERATURE  WHERE  PLOT  IS  TO  START 

C~  SCALE=MAX.  TEMPERATURE/40 

C  DD=LABEL  FOR  PLOT 

"SUBROUTINE  THEORY ( I Q , JQ, KQ , I  OP ) 
COMMON  GIVE,S,E,R,V,P,X,QSQ,NE,NP, IQ , JQ, KQ , I  OP , XTEST , NFQRM , NMQVE , A 

C0MM0N/C0MD£R/S4F,SLAT  " 

DIMENSION  GIVE( 300),S4F{ 120) 

OIMENS ION  DD( 10) 

DIMENSION  E(100),R{100),V{100),P(3,20),X<20),QSQ{3),DP(20),S(300) 

IF(NF0RM)1,1,2 
1  R2=0.19 

R3=7A.O 

ALPHA=i.020lE-2 
~~ALPHA2  =  ALPHA*  ALPHA" 

J  =  6 
"~SUM=J*{J+l)*(2*J+l)*T2»J-i)*(2*J*3>/5 

RM1=1.0    

~B=0.695056 

E2Q=0. 769854 


HCE=0.8610308E-5 

ANUM=-0.5*E2Q*R2«R3«ALPHA2*SUM»RMl/(  13.0»B) 
ALAT=0.5*E2Q*HCE*4.0 

RhAD(i>,  100)  TEMP,  SCALE 

"100  FORMAT (2F20. 5) 
2  ANM=P(KQ, 1 ) 
GAMMA=P(KQ,2 ) 
GO  TO  (3,4), IQ 
"3  ANM=P(JQ,1) 

GO  TO  5       

V"GAMMA=P(  JQ,  2  ) 

5  SLAT=-ALAT 
1F{ I0P)6, 7,8 

7  Ml=41 
GO  TO  9 

8  M1=NE  

GO  TO  9 

6  Ml  =  l 


9  DO  10  M=1,M1 

IF{ IOP) 11, 12,13 
11  T=V(NE) 

GO  TO  14 
"12~T  =  TEMP*SCAL£*FL0ATTMT~ 

GO  TO  14 

13  T=V(M) 

14  S4F(M) =-ANUM/T 
S(M)=-ANM»S4F(M)-GAMMA*SLAT 

10  CONTINUE 


-120- 

IF ( IOP)15, 16, 15 
16  P(2,2)=P(2,2)/P(2,i) 

WRITE(6,200)  P(2,l)  ,P(2,2) 
200  FGRMAT(26H0   SQRT(A20*»2  ♦  1/3A22»«2)=   F10.2,22H  "( 1-GAMMA) /( 1-SI 
1GMA)=   F10.2) 
READ(5, 101 )0D 

101  FURMATU0A6) 

WRITE (6,  102)00, TEMP, SCALE 

102  F0RMAT150H1   PLOT  OF  TEMPERATURE  VS  QUADRUPOLE  SPLITTING  IN  10A6/ 
'"140H0   NO.  OF  POINTS~~  =  41     IN  ITT  AL~TEMP~=  F10V1 ,~~~ 

218H   SCALE  FACTOR  »   F10.1) 
WRITE(6,103)  (S(  I  ),  1  =  1, Ml) 

103  F0RMAT(2X   20F6.3) 

CALL  J  PLOT ( S, A  1, A  1,0. 0,1 0.0, SCALE, TEMP, DO, IT" 

15  RETURN 

END  ~~ 


-121 


10 

20 
30 


THE  PARTIAL  DERIVATIVE  OF  THE  OUTPUT 
WITH  RESPECT  TO  THE  PARAMETER  IP  AT  THE 

CASE  OF  THE  SHIELDING  PARAMETERS 


DERI  V  2... 

THIS  FUNCTION  CALCULATES 

OF  THE  SUBROUTINE  THEORY 

POINT  JP. 

DERIV2  IS  FOR  THE  SPECIAL 

ONLY. 

FUNCTION  DERIVI IP, JP) 

DIMENSION  E1100),R(100),V(100),P(3,20),X(20),QSQ(3),S(300) 

DIMENSION  GIVE(300) 

DIMENSION  S4F1120) 

COMMON  GIVE,S,E,R,V,P,X,QSQ,NE,NP, 10 , JQ, KQ , I  OP, X  TEST , NFORM , MMOVE 

COMMON/COM0ER/S<tF,SLAT 

CALL  THEORY! 1,2,2, 1) 

IF(IP-l) 10,10,20 

DERIV  =  -S4F{  JPJ 

GO    TO    30 
DERIV=-SLAT 
RETURN 
END 
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